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Abstract:  The  U.S.  Army  Engineer  District,  Galveston,  has  recently 
enlarged  the  Houston  Ship  Channel  from  a  40-ft  (12.2-m)  depth  by  400-ft 
(122-m)  width  to  a  45-ft  (13.7-m)  depth  by  530-ft  (162-m)  width.  Previ¬ 
ously,  a  three-dimensional  (3-D)  numerical  model  study  was  implemented 
at  the  U.S.  Army  Engineer  Research  and  Development  Center’s  Coastal 
and  Hydraulics  Laboratory  to  evaluate  the  salinity  and  circulation  impact 
of  this  enlargement.  The  enlarged  channel  is  nearly  complete.  Prelimi¬ 
nary  evaluations  indicate  a  higher  than  anticipated  rate  of  deposition  in 
the  channel  reach  near  Atkinson  Island.  The  purpose  of  this  investigation 
was  to  determine  if  this  new  shoaling  rate  was  a  permanent  feature  of  this 
new  channel  or  a  single  anomalous  event.  If  this  considerably  larger  depo¬ 
sition  is  a  permanent  condition  then  the  planned  disposal  locations  must 
be  expanded.  The  project  will  determine  the  possible  effects  on  sedimen¬ 
tation  in  the  channel  based  on  dredging  records,  field  data,  vessel  effects, 
and  hydrodynamic  changes  due  to  the  channel’s  enlargement. 
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Preface 

This  report  represents  the  findings  of  an  investigation  of  the  possible 
causes  of  increased  shoaling  in  the  Houston  Ship  Channel  after  its  dimen¬ 
sions  were  increased  from  40  ft  x  400  ft  (12.2  m  x  122  m)  to  45  ft  x  530  ft 
(13.7  m  x  162  m).  This  unexpected  increase  in  shoaling  may  subside  in 
time  if  due  to  the  dredging  itself  or  the  channel  side  slopes  reestablishing 
equilibrium.  It,  however,  may  be  due  to  the  larger  channel  dimensions.  In 
which  case  the  increased  shoaling  would  be  a  permanent  condition.  Addi¬ 
tional  disposal  capability  would  then  be  required. 

This  investigation  was  conducted  from  August  2004  through  August  2005 
at  the  U.S.  Army  Engineer  Research  and  Development  Center  (ERDC)  by 
Dr.  R.  C.  Berger,  A.  R.  Carrillo,  L.  M.  Lee,  and  J.  N.  Tate  of  the  Coastal  and 
Hydraulics  Laboratory  (CHL).  The  data  collection  was  conducted  20- 
22  August  2004.  In  addition  to  Dr.  Berger,  the  data  collection  group  in¬ 
cluded  J.  R.  Bull,  H.  A.  Benson,  and  Dr.  K.  R.  Barry.  Funding  was  pro¬ 
vided  by  the  U.S.  Army  Engineer  District,  Galveston. 

The  work  was  performed  under  the  general  direction  of  Thomas  W. 
Richardson,  Director,  CHL,  Dr.  William  D.  Martin,  Deputy  Director,  CHL, 
B.  A.  Ebersole,  Chief,  Flood  and  Storm  Protection  Division,  CHL,  and  Dr. 
R.  T.  McAdory,  Chief,  Estuarine  Engineering  Branch,  CHL. 

COL  Richard  B.  Jenkins  was  Commander  and  Executive  Director. 

Dr.  James  R.  Houston  was  Director  of  ERDC. 
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1  Introduction 

Background 

The  U.S.  Army  Engineer  District,  Galveston,  recently  enlarged  the  Hous¬ 
ton  Ship  Channel  from  a  40-ft  (12.2-m)  depth  by  400-ft  (122-m)  width  to  a 
45-ft  (13.7-m)  depth  by  530-ft  (162-m)  width.  Previously,  a  three- 
dimensional  (3-D)  numerical  model  study  was  implemented  at  the  U.S. 
Army  Engineer  Research  and  Development  Center’s  (ERDC’s)  Coastal  and 
Hydraulics  Laboratory  (CHL)  to  evaluate  the  salinity  and  circulation  im¬ 
pact  of  this  enlargement.  In  Berger  et  al.  (1995a)  the  model  was  shown  to 
represent  the  salinity  and  circulation  in  the  earlier  channel  configuration. 
Then  in  Berger  et  al.  (1995b),  it  was  used  to  predict  the  impact  of  the 
enlarged  channel.  This  model  was  again  used  in  Carrillo  et  al.  (2002)  to 
evaluate  the  addition  of  barge  lanes  along  the  ship  channel  flanks.  The 
enlarged  channel  is  nearing  completion,  and  preliminary  evaluations  indi¬ 
cate  a  higher  than  anticipated  rate  of  deposition  in  the  channel  reach  near 
Atkinson  Island.  The  purpose  of  this  investigation  was  to  determine  if  this 
new  shoaling  rate  was  a  permanent  feature  of  this  enlarged  channel  or  a 
single  anomalous  event.  If  this  considerably  larger  deposition  is  a  perma¬ 
nent  condition,  then  the  planned  disposal  locations  must  be  expanded. 

The  historical  record  shows  that  once  the  bay  channel  became  deeper  than 
the  natural  depth  over  its  entire  length  (mid- 19 20’s),  the  dredging  rate  has 
been  fairly  constant.  This  may  mean  that  the  material  moves  into  the 
channel  from  the  sides  and  is  simply  moved  about  along  the  channel  by 
the  currents  and  vessel  entrainment,  or  this  could  be  the  result  of  other 
conditions.  The  historical  rate  of  dredging  within  the  bay  segment  of  the 
channel  has  been  about  2  million  cu  yd  (765,000  cu  m)  per  year.  The 
model  and  field  data  show  that  the  residual  near-bed  currents  are  rela¬ 
tively  large  over  most  of  the  southern  part  of  the  bay  channel,  approxi¬ 
mately  from  the  Texas  City  dike  up  to  about  where  the  Bayport  Channel 
intersects  the  main  channel.  The  residual  near-bed  currents  in  the  chan¬ 
nel  near  Atkinson  Island  are  small.  Over  the  entire  length  of  the  bay 
channel,  the  near-bed  currents  are  flood-directed.  That  is,  they  are  di¬ 
rected  landward  rather  than  toward  the  ocean.  Any  material  that  falls  into 
the  southern  part  of  the  channel  will  tend  to  move  upstream  to  the  region 
near  Atkinson  Island.  If  the  currents  are  not  strong  enough,  the  material 
may  deposit  before  it  reaches  this  area.  The  vessel  traffic  along  the  Hous- 
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ton  Ship  Channel  can  cause  resuspension  even  if  the  tidal  currents  are 
weak. 

The  preliminary  evaluation  of  the  required  dredging  quantity  for  the  new 
channel  configuration  appears  to  be  perhaps  4  million  cu  yd  (3  million 
cu  m)  per  year  in  just  the  channel  section  near  Atkinson  Island.  It  is  pos¬ 
sible  that  this  additional  dredging  quantity  is  material  that  has  escaped 
from  the  disposal  sites  or  low  density  fluff.  It  may  also  be  due  to  the  ad¬ 
justment  of  the  channel  side  slopes  to  the  new  configuration.  However,  it 
could  also  be  that  the  hydrodynamics  of  the  channel  and  bay  have  changed 
enough  to  result  in  this  change. 

Another  possible  reason  that  the  dredging  quantity  could  have  increased  is 
from  the  effect  of  vessel  traffic.  This  may  include  vessel-induced  erosion 
of  the  surrounding  shallows.  With  a  larger  channel,  the  vessels  can  move 
faster.  The  vessel  speed  is  limited  by  the  ability  of  the  water  to  move  out  of 
the  way  and  return  behind  the  vessel.  The  wave  speed  in  the  shallows 
around  the  vessel  may  be  slower  than  the  vessel.  In  this  circumstance,  the 
condition  is  like  that  of  supercritical  flow  in  open  channel  flow.  The  draw¬ 
down  wave  results  in  high  currents  and  substantial  shock  waves  moving 
along  and  behind  the  vessel.  It  is  possible  that  this  is  causing  more  sedi¬ 
ment  to  move  into  the  channel.  Also,  a  vessel  moving  in  the  channel  can 
reduce  the  deposition  rate  or  cause  erosion.  The  enlarged  channel  may 
have  reduced  the  erosion  potential  of  the  passing  vessels. 

Approach 

The  study  focuses  on  determining  what  the  ultimate  deposition  rate 
should  be  in  the  navigation  channel,  especially  in  the  Atkinson  Island 
reach,  for  the  45-  x  530-ft  (13.7-  x  162-m)  channel  (hereafter  referred  to  as 
the  "plan"  condition).  The  original  channel  configuration  of  40  x  400  ft 
(12.2  x  122  m)  is  termed  the  "base"  condition.  Several  different  compo¬ 
nents  of  the  study  are  undertaken. 

The  initial  step  is  to  analyze  the  dredging  records  obtained  from  the  Corps 
of  Engineers  District,  Galveston.  These  will  be  used  to  determine  the  pat¬ 
tern  of  deposition  before  the  enlargement  and  any  correlation  with  flow 
rate. 

Next,  field  data  collection  will  be  conducted: 


a.  To  provide  sediment  parameters  for  the  model. 
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b.  To  try  to  determine  the  source  of  the  extensive  sedimentation  in  the 
Atkinson  Island  reach  of  the  Houston  Ship  Channel. 

Sediment  data  include  bed  and  water  column  samples.  These  data  are 
used  to  estimate  critical  shear  stress,  settling  velocity,  grain-size  distribu¬ 
tion,  etc.  In  addition  to  characterizing  sediment  along  the  channel,  nearby 
shallows  are  also  sampled.  A  preliminary  assessment  of  resuspension  due 
to  ship  induced  currents  is  made  as  well.  This  is  used  to  determine  the 
potential  reduction  in  net  shoaling  in  the  base  and  enlarged  channel  con¬ 
ditions  due  to  the  presence  of  ships. 

The  next  step  is  to  analyze  the  currents  from  the  model  for  both  channel 
configurations  -  before  and  after  the  channel  enlargement.  The  average 
residual  currents  can  provide  an  indicator  of  the  shear  stress  vector  for  the 
before  and  after  enlargement  conditions.  The  questions  to  be  answered 
are  as  follows: 

a.  Is  there  substantially  larger  erosion  of  the  surrounding  shallows  with 
the  increased  channel  size? 

b.  Is  there  a  substantial  decrease  in  shear  stress  in  the  channel? 

c.  Do  the  residual  currents  suggest  a  large  shift  in  the  distribution  of 
sediment  similar  to  that  noticed  for  this  first  maintenance  cycle? 

The  final  step  is  to  improve  the  previous  3-D  model  to  include  more  de¬ 
tailed  representation  of  areas  that  are  of  interest  for  this  particular  study. 
The  model  is  then  used  to  analyze  the  impact  of  the  channel  enlargement 
by  introducing  a  sediment  tracer  in  the  model  channel  for  the  40-  x  400-ft 
(12.2-  x  122-m)  channel  and  the  “as  built”  45-  x  530-ft  (13.7-  x  162-m)  con¬ 
ditions.  This  gives  an  idea  of  the  change  one  should  see  in  the  dredging 
distribution  and  an  indication  of  any  change  in  trapping  efficiency.  This 
reflects  only  the  impact  of  the  change  in  hydrodynamics  and  therefore  is 
only  a  rough  estimate  of  the  expected  permanent  change  in  shoaling. 

The  analysis  of  the  dredging  data  is  in  the  second  chapter.  The  third  chap¬ 
ter  discusses  the  data  collection,  the  fourth  chapter  discusses  the  vessel- 
induced  effects  on  sedimentation,  and  the  fifth  chapter  details  model  re¬ 
sults.  First  the  impact  of  the  hydrodynamic  change  due  to  the  channel 
enlargement  and  its  likely  effect  on  sedimentation  is  investigated.  Second, 
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the  model  is  run  with  a  sediment  tracer.  The  final  chapter  is  a  summary 
with  conclusions  from  the  entire  report. 
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2  Dredge  Data  Analysis 

The  dredging  records  of  the  Houston  Ship  Channel  for  the  period  in  which 
the  channel  was  in  the  40-ft  (12.2-m)  deep  by  400-ft  (122-m)  wide  condi¬ 
tion  are  analyzed.  There  are  two  goals  of  this  analysis:  a)  Determine  the 
distribution  of  sedimentation  during  the  period  in  which  the  channel  was 
in  the  40-  x  400-ft  (12.2-  x  122-m)  configuration.  This  is  to  be  used  as 
validation  data  for  the  sediment  model;  which  will  be  reported  on  later, 
b)  Determine  the  change  in  shoaling  over  time  during  this  same  period, 
and  attempt  to  find  a  causal  relationship  with  freshwater  inflow  to  the  sys¬ 
tem  or  other  factors. 

Figure  1  shows  the  Houston  Ship  Channel  area  and  locations  of  interest  for 
this  analysis. 


Figure  1.  Houston  Ship  Channel  area  map. 


Historical  dredging  records  were  received  from  the  U.S.  Army  Corps  of 
Engineer  District,  Galveston,  in  EXCEL  spreadsheet  form  for  all  reaches 
along  the  Houston  Ship  Channel  over  a  time  period  from  1949  to  2003. 
These  data  were  sorted,  and  the  reaches  of  Bolivar  Roads  to  Red  Fish  Reef 
and  Red  Fish  Reef  to  Morgan’s  Point  were  separated  from  the  other 
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reaches.  The  data  were  provided  as  total  yardage  dredged  on  a  given  date 
for  a  given  station  range.  The  stations  were  given  as  an  upstream  station 
number  and  a  downstream  station  number  for  each  yardage  total,  with 
station  numbers  being  1000  ft  (  302.8  m)  increments.  The  station  num¬ 
bers  were  used  as  the  guide  to  join  the  two  reaches  together  into  one  data 
set.  Bolivar  Roads  was  taken  as  sta  130,  with  the  station  numbers  decreas¬ 
ing  as  one  travels  upstream  to  Morgan’s  Point  and  sta  o  (see  Figure  2). 

The  analysis  presented  here  is  to  determine  if  a  correlation  can  be  made 
between  the  dredging  records  and  the  freshwater  flows  into  the  system 
during  this  time. 


Figure  2.  Houston  Ship  Channel  station  map  (emergent  features  are  shaded  green). 


In  order  to  determine  a  yardage  for  each  station,  the  number  of  stations 
included  in  each  dredging  segment  was  divided  into  the  total  yardage  for 
that  segment.  The  yardage  was  summed  over  time  for  each  of  the  stations 
and  averaged  by  the  total  amount  of  time,  -34  years.  Note  that  the  time  of 
interest  for  this  project  is  after  the  navigation  channel  was  dredged  to  40  ft 
(12.2  m).  It  was  determined  that  this  deepening  was  performed  by  Fiscal 
Year  1964.  The  dredge  data  prior  to  this  deepening  were  not  included  in 
the  analysis  presented  here.  These  sections  were  last  dredged,  for  mainte¬ 
nance  only,  in  1998,  to  maintain  the  enlarged  -  widened  and  deepened  - 
channel.  This  manipulation  allowed  for  the  determination  of  an  average 
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shoaling  rate  per  station  over  the  entire  dredging  period  for  each  of  the 
stations.  The  average  shoaling  rate  per  station  is  given  in  Figure  3. 


Average  Shoaling  Rate  Since  40  ft  deepening 

total  yardage  summed  since  1963  dredging  divided  by  total  time  from  1963  dredge  to  1998  dredge 
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□  Average  Shoaling  Rate  Since  40  ft  deepening 


Figure  3.  Average  shoaling  rate  per  station  (total  cubic  yards  per  total  number  of  years  for 

each  station). 


The  dredge  data  were  separated  by  reach  with  station  72  as  the  dividing 
line  for  Red  Fish  Reef.  The  total  yardage  for  each  dredge  occurrence  was 
summed  over  the  reach  (Morgan’s  Point  to  Red  Fish  Reef:  stations  0-72 
and  Red  Fish  Reef  to  Bolivar  Roads:  stations  72-130).  Each  yardage  sum 
was  averaged  by  the  number  of  stations  in  the  reach  and  divided  by  the 
time  since  the  previous  dredge  event  for  that  reach  in  order  to  obtain  an 
average  shoaling  per  year  over  the  reach.  The  average  shoaling  rate  per 
reach  over  time  is  shown  in  Figure  4  for  Morgan’s  Point  to  Red  Fish  Reef 
and  in  Figure  5  for  Red  Fish  Reef  to  Bolivar  Roads. 

Results  of  this  analysis  showed  that  the  shoaling  rate  was  higher  in  the 
upper  half  of  the  bay  channel  for  the  40-  x  400-ft  (12.2-  x  122-m)  condition 
(this  is  in  the  region  north  of  Redfish  Reef  and  along  Atkinson  Island).  In 
this  upper  reach  the  sedimentation  dredging  volume  dropped  dramatically 
and  fairly  consistently  over  the  34  years  of  this  channel  configuration.  The 
dredging  rate  in  the  recent  years  was  only  about  one-fourth  of  that  re¬ 
quired  immediately  after  the  enlargement  to  40  x  400  ft  (12.2  x  122  m). 
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The  decline  is  steady  throughout  the  34  years.  In  the  lower  bay  there  was 
no  strong  trend  over  this  period. 


Average  Shoaling  Rate  per  Year  (Morgans's  Point  to  Red  Fish  Reef  -  station  0  to  station  72) 
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Figure  4.  Average  shoaling  rate  for  Morgan's  Point  to  Red  Fish  Reef  (total  cubic  yards  per  station  in  the 

reach  per  year  between  dredging). 
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Average  Shoaling  Rate  per  Year  (Red  Fish  Reef  to  Bolivar  Roads  -  station  72  to  station130) 


Years 

■  Avg  Shoaling  Rate  over  reach 

Figure  5.  Average  shoaling  rate  for  Red  Fish  Reef  to  Bolivar  Roads  (total  cubic  yards  per  station  in  the 

reach  per  year  between  dredging). 

Stream  flow  data  for  the  Trinity  River  and  the  San  Jacinto  River  over  the 
same  time  period  as  the  dredging  records  were  obtained  from  the  U.S. 
Geological  Survey  (USGS)  Web  site.  In  order  to  get  data  over  the  entire 
period  needed,  data  sets  were  taken  from  gauges  upstream  of  the  rivers’ 
entrances  to  the  Houston  Ship  Channel  model  domain.  Trinity  River  data 
were  taken  from  the  gauge  at  Liberty,  TX,  while  San  Jacinto  data  were 
taken  from  the  western  fork  at  Conroe,  TX,  and  the  eastern  fork  at  Cleve¬ 
land,  TX.  The  purpose  of  the  stream  flow  data  was  to  determine  if  a  corre¬ 
lation  exists  between  it  and  the  dredging  records;  therefore  a  general  trend 
was  only  necessary  at  this  point  and  data  from  gauges  slightly  upstream 
were  reasonable.  The  flow  data  are  presented  in  Figures  6-8. 
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San  Jacinto  Discharge  @  Conroe,  Tx  (1965-1998) 
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Figure  6.  Flow  data  for  San  Jacinto  River  western  fork. 
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It  was  thought  that  perhaps  the  freshwater  inflow  decreased  dramatically 
leading  to  this  drop  in  sedimentation.  However,  the  flow  data  do  not  sup¬ 
port  this  hypothesis  since  the  data  showed  somewhat  more  average  flow 
during  the  years  of  less  dredging. 

It  is  possible  that  the  channel  stabilized  over  time  causing  less  sediment  to 
erode  from  its  banks.  Perhaps  dredging  practices  improved  so  that  less 
dredged  material  was  allowed  to  reenter  the  channel.  Another  possibility 
for  the  drop  in  dredging  over  this  period  is  the  change  in  the  sediment  load 
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San  Jacinto  Discharge  @  Cleveland,  Tx  (1965-1998) 


Time  (days) 


Figure  7.  Flow  data  for  San  Jacinto  River  eastern  fork. 


Trinity  Discharge  @  Liberty,  Tx  (1965-1998) 


Figure  8.  Flow  data  for  Trinity  River. 
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from  upstream  sources  to  the  bay.  Ward  and  Armstrong  (1993)  found  that 
the  suspended  solids  in  the  upper  Houston  Ship  Channel  had  dropped  to 
one-third  of  that  from  25  years  earlier.  They  attribute  this  reduction  to 
improved  waste  treatment,  altered  land  use,  and  impoundments  on  the 
principal  rivers.  Solis  et.  al.  (1994)  found  that  the  Trinity  River  gauge  at 
Romayor,  TX  (USGS  gauge  08066500  located  approximately  74  km  (46 
miles)  upstream  from  the  river’s  entrance  to  the  bay)  shows  a  significant 
reduction  in  downstream  sediment  inputs  after  the  completion  of  Lake 
Livingston  in  1968.  White  et.  al.  (2002)  indirectly  confirmed  this  finding 
by  noting  that  the  Trinity  River  delta  was  prograding  through  most  of  the 
20th  century,  with  a  transition  to  degradation  beginning  sometime  be¬ 
tween  1956  and  1974. 

There  actually  is  a  decline  in  total  suspended  sediment  throughout  most  of 
the  bay.  The  goal  of  this  analysis  was  to  study  the  link  between  dredging 
quantities  and  the  freshwater  inflow  to  the  system.  There  does  not  appear 
to  be  a  clear  link.  It  appears  that  the  perception  that  the  historical  shoal¬ 
ing  rate  in  the  channel  had  not  increased  with  channel  enlargement  may  in 
fact  be  due  to  reduced  sediment  load  caused  by  other  factors,  such  as  the 
impoundment  of  the  rivers,  waste  treatment,  altered  land  use,  or  improved 
dredge  disposal  practices. 
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3  Field  Data  Collection  and  Analysis 

The  ERDC-CHL  Measurement  and  Analysis  Group  sent  a  team  to  the 
Houston  Ship  Channel  -  Galveston  Bay  area  to  gather  field  data  from  20 
to  22  August  2004.  A  summary  of  observations  is  given  in  this  chapter 
along  with  the  actual  data  and  analysis  performed. 

Generally  the  channel  was  found  to  contain  significant  deposits  along  the 
edges  with  the  center  maintaining  the  49  ft.  (14.9-m)  or  so  depth  dredged 
during  channel  construction.  When  not  meeting  or  passing  other  vessels 
the  ship’s  path  was  along  the  centerline.  It  appears  that  the  forces  (return 
currents  and  propeller)  are  enough  to  keep  deposits  away  from  the  chan¬ 
nel  centerline. 

Bed  samples  were  collected  using  a  3-in.-  (7.62-cm-)  diam  drag  bucket  bed 
sampler  and  1-1/2-in.-  (3.81-cm-)  diam  push  core  samples.  Many  more 
samples  were  collected  than  were  analyzed.  In  the  case  of  the  push  core 
samples  many  were  collected  to  note  the  type  of  material  present.  Analysis 
of  a  few  representative  samples  was  then  made,  rather  than  the  entire 
group.  Suspended  sediment  concentrations  were  collected  using  a  3-1/2- 
in.-  (8.89-cm-)  diam  Niskin  tube.  The  complete  analysis  data  for  calculat¬ 
ing  suspended  sediment  concentrations,  fall  velocity,  erodibility  parame¬ 
ters,  bulk  density,  and  sediment  gradations  are  supplied  in  the  appendix. 
This  chapter  summarizes  the  findings. 

The  shallows  on  the  east  side  of  the  channel  in  the  Atkinson  Island  reach 
were  of  coarser  material  than  on  the  west  side  at  about  the  same  distance 
from  the  channel  and  depth.  This  can  be  seen  from  the  data  in  Table  1  and 
Figure  9  that  give  the  fraction  of  the  sample  that  is  silt/clay  and  that  which 
is  coarser  (here  it  is  referred  to  as  sand).  The  samples  at  090009, 112445, 
and  123924  on  the  east  side  of  the  channel  show  silt/clay  fractions  of  less 
than  10  percent.  While  on  the  west  side  of  the  channel,  stations  114554, 
121418,  and  123028  have  silt/clay  fractions  of  around  20  to  35  percent. 

The  east  side  of  the  channel  is  exposed  to  somewhat  greater  wave  forces 
due  to  the  reflection  caused  by  the  presence  of  Atkinson  Island.  This 
would  be  true  for  vessel-produced  waves  in  particular. 

On  the  west  side  of  the  channel  the  bed  material  becomes  progressively 
finer  as  one  moves  away  from  the  channel.  The  progression  from  stations 
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121418,  082124,  and  074151  results  in  silt/clay  fractions  of  about  15,  55, 
and  90  percent,  respectively.  This  could  be  the  usual  delta  formation  of  a 
river  moving  into  a  bay,  the  larger  sizes  settling  out  sooner  than  the  finer 
size  particles.  Of  course,  near  the  channel  the  vessel  return  currents  and 
waves  also  play  a  part. 


Table  1.  Sand-silt/clay  separation  analysis  results  giving  percentage  of 
each  at  several  points. 


Sand-Silt/Clay  Separation  Analysis 

PROJECT:  Houston  Ship  Channel 

Analysis  complete! 

9/22/04 

P.I.:  C.  Berger 

Analyst:  D.  B. 

Sampling  date:  8/20/04 

Date 

Range 

Time 

Depth 

Sed. 

Tare 

Gross 

Pan 

Net 

Total  sed. 

Sed. 

Sta.  # 

CST 

(ft.) 

type 

wt.  g 

wt.  g 

wt.  g 

wt.  g 

wt.  g 

% 

8/20/04 

090009 

10.0 

sand 

72.75 

83.10 

0.1991 

10.15 

11.21 

90.55 

silt/clay 

152.25 

153.11 

1.06 

9.45 

8/20/04 

093838 

48.0 

sand 

72.42 

72.50 

0.0163 

0.06 

3.20 

1.99 

silt/clay 

152.67 

155.79 

3.14 

98.01 

8/20/04 

112445 

6.0 

sand 

72.07 

80.79 

0.1743 

8.55 

9.44 

90.53 

silt/clay 

161.82 

162.54 

0.89 

9.47 

8/20/04 

114554 

11.0 

sand 

73.06 

78.99 

0.3016 

5.63 

7.29 

77.21 

silt/clay 

160.96 

162.32 

1.66 

22.79 

8/20/04 

121418 

13.0 

sand 

71.96 

81.39 

0.4552 

8.97 

10.75 

83.49 

silt/clay 

154.82 

156.14 

1.78 

16.51 

8/20/04 

123028 

11.0 

sand 

72.60 

79.85 

0.8610 

6.39 

9.95 

64.21 

silt/clay 

164.28 

166.98 

3.56 

35.79 

8/20/04 

123924 

9.0 

sand 

72.18 

80.33 

0.1235 

8.03 

8.62 

93.11 

silt/clay 

163.78 

164.25 

0.59 

6.89 

8/22/04 

Push  core  #1 

074151 

10.5 

sand 

88.31 

89.04 

0.1458 

0.58 

5.33 

10.96 

10/26/04 

silt/clay 

164.40 

169.00 

4.75 

89.04 

8/22/04 

Push  core  #3 

082124 

9.0 

sand 

88.31 

94.18 

1.0701 

4.80 

11.03 

43.52 

11/2/04 

silt/clay 

164.40 

169.56 

6.23 

56.48 

Now  consider  the  channel  bed  material.  Figures  10-14  show  stations 
093838, 110820, 114205, 120447,  and  123417,  respectively.  These  show 
the  sediment  gradations  by  volume  using  a  Coulter  counter.  The  order  of 
the  figures  begins  upland  in  the  bayou  section  of  the  Houston  Ship  Chan¬ 
nel  and  then  moves  south  along  the  channel  to  near  the  midreach  of  Atkin¬ 
son  Island.  These  plots  show  modes  at  about  100  microns  (0.0039  in.) 
(very  fine  sand),  20  microns  (0.00079  in.)  (silt),  and  approximately  4  mi¬ 
crons  (0.00016  in.)  (clays).  The  trend  moving  down  the  channel  through 
these  stations  is  a  shift  toward  coarser  material.  That  is,  the  upland  sta¬ 
tions  are  predominantly  clays  and  silts.  The  lower  stations  are  more  sand. 
This  trend  would  be  consistent  with  material  moving  upland  in  the  chan¬ 
nel  and  it  is  known  that  the  channel’s  net  drift  is  upland.  This  has  been 
noted  in  Bobb  et  al.  (1973);  Ward  (1980);  and  Berger  et  al.  (1995a).  These 
references  confirmed  that  over  most  of  the  lower  bay  channel,  the  drift  is 
in  the  flood  direction  (upland)  over  the  water  column.  From  Bobb  et.  al. 
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(1973),  at  about  the  southern  tip  of  Atkinson  Island,  the  channel  becomes 
more  stratified.  In  this  region  of  the  channel  and  upland,  the  bottom  cur¬ 
rents  are  usually  flood-directed  while  the  surface  currents  tend  to  be  ebb- 
directed.  The  radioactive  sediment  tracer  study  by  Hart  (1969)  demon¬ 
strated  that  the  sediments  were  moving  strongly  upland  (flood)  in  the 
channel. 


Figure  9.  Grain  class  separation  percentages  for  several  field  data  points 


The  open  bay  regions  in  the  vicinity  of  Atkinson  Island  (west  of  the  naviga¬ 
tion  channel  and  also  the  region  east  of  the  island)  contain  a  great  deal  of 
silts  and  clays.  The  areas  exposed  to  ship  traffic  are  sandier,  as  a  likely  re¬ 
sult  of  natural  armoring.  This  is  noticeable  at  the  northern  tip  of  Atkinson 
Island.  Along  the  shore  to  the  west  are  shells  and  sands.  To  the  east  is 
finer  material.  The  disposal  site  just  south  of  Atkinson  Island  is  more  firm 
material. 

Some  vessels  passing  along  the  waterway  created  bores  of  about  a  foot  in 
height  in  shallow  areas.  These  include  the  area  along  Hogg  Island  and  pe- 
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riodically  along  Atkinson  Island.  These  locations  were  shallow  and  the 
bed  material  was  found  to  be  firm.  Easily  erodable  material  would  have 
been  swept  away  by  the  wave. 

Stations  164015  and  163007  are  shown  in  Figures  15  and  16,  respectively. 
These  stations  are  in  the  channel  and  the  nearby  shallows.  Both  show  sig¬ 
nificant  fine  material.  There  is  a  midbay  marsh  site  being  created  nearby. 
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Coulter®  LS  Particle  Size  Analyzer 


File  name: 

093838 

Group  ID: 

Houston  Ship  Channel 

Sample  ID: 

093838  48  ft.  8-20-04 

Location: 

Range  1 ,  channel 

Operator: 

DB 

Comments: 

oxi.  3  times 

Optical  model: 

Fraunhofer,  rfz 

LS  100Q 

Fluid  Module 

Volume  Statistics  (Arithmetic)  093838 
Calculations  from  0.375  urn  to  948.3  pm 


Volume: 

100% 

Mean: 

9.679  pm 

S.D.: 

8.424  pm 

Median: 

6.384  pm 

C.V.: 

87.0% 

D(3,2): 

2.852  pm 

Skewness 

0.658  Right  Skewed 

Mode: 

21.69  pm 

Kurtosis: 

-0.947  Pltykurtic 

%< 
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90 
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17.28 

22.67 

Figure  10.  Coulter  counter  analysis  for  Point  093838. 


ERDC/CHL  TR-06-8 


18 


Coulter®  LS  Particle  Size  Analyzer 


File  name: 

110820 

Group  ID: 

Houston  Ship  Channel 

Sample  ID: 

110820  39  ft.  8-20-04 

Location: 

Range  3,  channel,  Morgan’s  Point 

Operator: 

DB 

Comments: 

oxi.  3  times 

Optical  model: 

Fraunhofer,  rfz 

LS  100Q 

Fluid  Module 

Volume  Statistics  (Arithmetic)  110820 


Calculations  from  0.375  jam  to  948.3  pm 


Volume: 

100% 

Mean: 

1 1.30  pm 

S.D.: 

12.44  pm 

Median: 

5.914  pm 

C.V.: 

110% 

D(3,2): 

2.776  pm 

Skewness: 

1 .445  Right  Skewed 

Mode: 

4.047  pm 

Kurtosis: 

1  251  Leptokurtic 

%< 

10 

25 

50 

75 

90 

(uim 

0.990 

2.164 

5.914 

16.43 

31.62 

Figure  11.  Coulter  counter  analysis  for  Point  110820. 
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Coulter®  LS  Particle  Size  Analyzer 


File  name: 

114205 

Group  ID: 

Houston  Ship  Channel 

Sample  ID: 

114205  47  ft.  8-20-04 

Location: 

Range  4,  channel,  near  Mooring  Basin 

Operator: 

DB 

Comments: 

oxi.  3  times 

Optical  model: 

Fraunhofer,  rfz 

LS100Q 

Fluid  Module 

Differential  Volume  (Average)  (2  S.D.) 


Volume  Statistics  (Arithmetic) 

114205 

Range  4  L 

Channel  ■«*.  \ 

/  /X 

•  LO 

Near  mooring  \ 

.  i 

Calculations  from  0.375  pm  to  948.3  pm 

basin  ■. 

i 

\ 

Volume: 

100% 

/  A 

\ 

Mean: 

8.410  pm 

S.D.: 

8.922  pm 

*  V 

Median: 

4.714  pm 

C.V.: 

106% 

/ 

D(3,2): 

2.503  pm 

Skewness: 

1.471  Right  Skewed 

A  \ 

Mode: 

4.877  pm 

Kurtosis: 

1.417  Leptokurtic 

V 

\  \ 

—  \ 

\ 

■ 

* 

%<  10 

25  50 

75  90 

/ 

j 

■ 

pm  0.925 

1.940  4.714 

11.80  22.91 

/ 

■ 

/ 

** 

Figure  12.  Coulter  counter  analysis  for  Point  114205. 
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Coulter®  LS  Particle  Size  Analyzer 


File  name: 

120447 

Group  ID: 

Houston  Ship  Channel 

Sample  ID: 

120447  37  ft.  8-20-04 

Location: 

Range  5,  channel 

Operator: 

DB 

Comments: 

oxi.  3  times 

Optical  model: 

Fraunhofer,  rfe 

LS  100Q 

Fluid  Module 

Differential  Volume  (Average)  (2  S.D.) 


Particle  Diameter  (pm) 


Volume  Statistics  (Arithmetic)  120447 
Calculations  from  0.375  (im  to  948.3  pm 


Volume: 

Mean: 

Median: 

D(3,2): 

Mode: 

%<  10 

pm  1 .386 


100% 

68.61  pm 
55.97  pm 
4.623  pm 
116.3  pm 

25  50 

4.311  55.97 


S.D.: 

C.V.: 

Skewness: 

Kurtosis: 

75  90 

115.3  147.8 


67.74  pm 
98.7% 

0.899  Right  Skewed 
0.429  Leptokurtic 


Figure  13.  Coulter  counter  analysis  for  Point  120447. 
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Coulter®  LS  Particle  Size  Analyzer 


File  name: 

123417 

Group  ID: 

Houston  Ship  Channel 

Sample  ID: 

123417  39  ft.  8-20-04 

Location: 

Range  6,  channel 

Operator: 

DB 

Comments: 

oxi.  3  times 

Optical  model: 

Fraunhofer  rfz 

LS100Q 

Fluid  Module 

Differential  Volume  (Average)  (2  S.D.) 


Particle  Diameter  (pm) 


Volume  Statistics  (Arithmetic)  123417 
Calculations  from  0.375  jam  to  948.3  pm 


Volume: 

Mean: 

Median: 

D(3,2): 

Mode: 


100% 
57.19  pm 
55.14  pm 
5.072  pm 
96.49  jim 


S.D.: 

C.V.: 

Skewness: 

Kurtosis: 


50.91  pm 
89% 

0.483  Right  Skewed 
-0.801  Platykurtic 


%<  10  25  50  75  90 

pm  1.567  5.135  56.14  96.51  126.2 


Figure  14.  Coulter  counter  analysis  for  Point  123417. 
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Coulter®  LS  Particle  Size  Analyzer 


File  name: 

164015 

Group  ID: 

Houston  Ship  Channel 

Sample  ID' 

164015  51  ft.  8-21-04 

Location: 

Channel,  near  marsh  site 

Operator: 

DB 

Comments: 

oxi.  3  times 

Optical  model: 

Fraunhofer,  rfz 

LS  100Q 

Fluid  Module 

Differential  Volume  (Average)  (2  S.D.) 


Particle  Diameter  (pm) 


Volume  Statistics  (Arithmetic)  164015 
Calculations  from  0.375  pm  to  948.3  pm 


Volume: 

100% 

Mean: 

11.11  pm 

S.D.: 

14.18  pm 

Median: 

5.314  pm 

C.V.: 

128% 

D(3,2): 

2.657  pm 

Skewness 

2.460  Right  Skewed 

Mode: 

4.443  pm 

Kurtosis: 

7.947  Leptokurtic 

%< 

10 

25 

50 

75 

90 

pm 

0.959 

2.054 

5.314 

14.88 

29.56 

Figure  15.  Coulter  counter  for  Point  164015. 
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Coulter®  LS  Particle  Size  Analyzer 


File  name 

163007 

Group  ID: 

Houston  Ship  Channel 

Sample  ID: 

163007  13  ft.  8-21-04 

Location: 

East  shallows,  near  marsh  site 

Operator: 

DB 

Comments: 

oxi.  3  times 

Optical  model: 

Fraunhofer,  rfz 

LS  100Q 

Fluid  Module 

Differential  Volume  (Average)  (2  S.D.) 


Particle  Diameter  (pm) 


Volume  Statistics  (Arithmetic)  163007 


Calculations  from  0.375  pm  to  948.3  |nm 


Volume: 

100% 

Mean: 

14.15  jam 

S.D.: 

22.42  pm 

Median: 

5.654  pm 

C.V.: 

158% 

D(3,2): 

2.784  pm 

Skewness: 

3.072  Right  Skewed 

Mode: 

4.443  pm 

Kurtosis: 

10.23  Leptokurtic 

%< 

10 

25 

50 

75 

90 

|j,m 

1.001 

2.182 

5.654 

15.97 

34.14 

Figure  16.  Coulter  counter  for  Point  163007. 
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Surveys  of  the  region  have  shown  that  material  spreads  a  considerable  dis¬ 
tance  from  the  midbay  marsh  site  and  appears  to  be  entering  the  channel, 
so  these  fines  may  be  from  that  location.  However,  they  could  be  from  the 
large  sources  to  the  east  and  west  of  the  channel,  and  unrelated  to  this 
marsh  site. 
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4  Vessel  Effects  on  Sediment  Suspension 


In  this  chapter,  the  influence  of  the  enlarged  channel  on  the  speeds  attain¬ 
able  by  vessels  in  the  Houston-Galveston  Navigation  Channel  is  examined. 
This  investigation  compares  the  same  vessels  in  the  previous  and  current 
configurations  using  a  typical  cross  section  for  each.  The  year  1990  con¬ 
figuration  of  the  bay  part  of  the  Houston  Ship  Channel  consists  of  a  40-ft 
(12.2-m)  mlw  depth  and  a  400-ft  (122-m)  base  width.  It  is  assumed  that 
the  side  slopes  are  1  Vertical  ^Horizontal  up  to  a  depth  at  mlw  of  10  ft 
(3  m).  For  the  enlarged  channel  (which  is  the  present  configuration)  the 
depth  is  45-ft  (13.7-m)  mlw  and  530-ft  (162-m)  base  width.  The  same  side 
slope  and  ultimate  depth  at  the  shoulder  of  10-ft  (3-m)  mlw  is  used.  The 
first  section  estimates  the  drawdown  and  return  currents  due  to  moving 
vessels,  then  describes  this  impact  to  the  waterway.  The  second  section 
estimates  the  erosion  due  to  vessel  passage. 

Vessel  Hydrodynamics  in  Houston  Ship  Channel 

The  equations  developed  by  Jansen  and  Shijf  (1953)  are  used  to  calculate 
the  cross-section  average  return  current  and  drawdown.  The  section  con¬ 
taining  just  the  channel  center  section  and  the  side  slopes  is  treated  as  the 
cross  section  and  does  not  consider  the  surrounding  shallows  of  the  bay. 
This  equation  is  actually  just  Bernoulli’s  energy  equation  about  the  vessel 
along  with  conservation  of  mass.  The  two  equations  are  as  follows: 


z 


(V+VJ2-V- 


2  9 


VsA={Vr+Vs)Aw 


(1) 


where: 

z  =  the  drawdown 

Vs  =  the  speed  of  the  vessel  relative  to  a  fixed  reference  frame 

Vr  =  the  return  current 

g  =  the  acceleration  due  to  gravity 

Ac  =  the  cross-sectional  area  in  front  of  the  vessel 
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Aw  =  the  cross-sectional  area  at  midlength  of  the  vessel. 

This  method  assumes  uniformity  over  a  cross-section,  no  friction  losses, 
and  zero  ambient  currents. 

Maynord  (2005b)  recorded  vessel  speeds  along  with  water  data  in  the 
channel  near  Atkinson  Island.  This  vessel  data  will  be  used  as  an  indicator 
of  vessels  that  demonstrated  strong  hydraulic  impacts.  Table  2  contains 
information  on  these  vessels. 


Table  2.  Description  of  vessel  dimensions  found  to  be  hydraulically  significant  in  Houston 

Ship  Channel  near  Atkinson  Island. 


Boat  Number 

Length,  ft  (m) 

Beam,  ft  (m) 

Draft,  ft  (m) 

1 

760  (232) 

106(32.3) 

30.5  (9.30) 

2 

760  (232) 

140  (42.7) 

30.5  (9.30) 

3 

899  (274) 

180  (54.9) 

34.6  (10.5) 

This  section  calculates  the  maximum  ship  velocity,  Vmax ,  for  the  two  chan¬ 
nels.  This  is  the  maximum  speed  attainable,  however,  this  is  fuel  ineffi¬ 
cient.  The  speed  that  is  typically  used  is  about  85  percent  of  the  maxi¬ 
mum.  Table  3  records  the  drawdown,  return  current,  and  the  depth  hcr  at 

which  the  Froude  number  is  1  for  this  ship  speed.  hcr  then  is  just  equal  to 

YLm 

g  ' 


Table  3.  Effect  of  channel  size  on  vessel  speed,  drawdown,  return  current,  and  critical  depth. 


Boat 

Channel,  ft  (m) 

y max  A/S 

(m/s) 

85  percent  Vmax 

Vs,f/s 

(knots) 

z ,  ft  (m) 

Vr,f/s 

(m/s) 

hcr  ’ ft  (m) 

1 

40X400 

(12.2X122) 

17.3  (5.27) 

14.7 

(8.7) 

1.5  (0.46) 

3.0  (0.91) 

6.7  (2.04) 

45X530 

(13.7X162) 

20.4(6.22) 

17.3  (10.2) 

1.4 

(0.43) 

2.5  (0.76) 

9.3  (2.84) 

2 

40X400 

(12.2X122) 

15.4  (4.70) 

13.1 

(7.8) 

1.6 

(0.49) 

3.5(1.07) 

5.3  (1.62) 

45X530 

(13.7X162) 

18.5  (5.64) 

15.7 

(9.3) 

1.6 

(0.49) 

2.9  (0.88) 

7.7  (2.35) 

3 

40X400 

(12.2X122) 

12.4  (3.78) 

10.5 

(6.2) 

1.6 

(0.49) 

4.2  (1.28) 

4.8  (1.46) 

45X530 

(13.7X162) 

15.5  (4.73) 

13.2 

(7.8) 

1.7 

(0.52) 

3.6  (1.10) 

5.4  (1.65) 
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The  important  point  of  this  table  is  that  the  larger  channel  allows  the  ves¬ 
sels  to  travel  faster.  The  vessel  can  travel  faster  in  the  channel  than  a  free 
surface  wave  can  travel  in  the  surrounding  shallows.  If  the  shallows  have  a 
depth  less  than  hcr  then  a  bore  will  be  formed.  The  enlarged  channel  low¬ 
ers  the  return  current  magnitude  and  had  little  impact  on  the  drawdown. 
However,  since  the  vessels  can  move  faster  in  the  larger  channel,  the  criti¬ 
cal  depth  is  much  deeper.  Therefore  ships  will  produce  bores  from  the 
drawdown  over  a  number  of  the  surrounding  shallows. 

The  importance  of  this  effect  is  that  the  measured  shoaling  has  shown  a 
significant  increase  since  the  channel  was  deepened.  Since  the  deepened 
channel  has  only  been  in  place  a  short  time,  there  is  little  certainty 
whether  this  is  going  to  be  the  typical  shoaling  rate  or  if  this  is  within  the 
variability  found  in  dredging  requirements.  There  is  ample  evidence  that 
there  is  erosion  of  shoreline  and  shallow  regions,  which  seems  to  have 
been  going  on  even  before  the  channel  was  enlarged.  Examples  are  the 
erosion  of  Goat  Island  (Maynord  2000)  and  also  the  islands  separating  the 
barge  mooring  basin  (Maynord  2005a)  from  the  main  channel.  If  the  in¬ 
creased  deposition  results  from  vessel-induced  shoreline  erosion,  then  one 
would  expect  this  to  decline  over  time  as  the  nearby  areas  are  eroded  and 
naturally  armored. 

Effects  of  Vessels  on  Shoaling 

As  a  vessel  moves  along  the  waterway,  it  creates  a  drawdown  and  a  return 
current  around  the  vessel.  The  return  currents  will  create  a  stress  on  the 
bed  that  can  cause  erosion.  These  periods  of  erosion  will  reduce  the  longer 
term  channel  deposition.  The  increased  channel  size  could  decrease  the 
stress  caused  by  the  vessel  passage  and  result  in  an  apparent  increase  in 
channel  shoaling.  The  propeller  may  also  create  erosion.  However,  in  this 
analysis  only  cross-sectional  effects  of  the  vessel  body  are  determined  (not 
including  the  propeller  effects).  To  perform  this  analysis,  the  estimated 
return  currents  for  the  largest  of  the  three  previously  mentioned  vessels, 
information  from  the  analysis  of  the  field  collected  sediment,  and  a  crude 
estimate  of  ship  traffic  are  utilized.  This  analysis  will  show  how  signifi¬ 
cant  the  vessel  traffic  may  be  on  the  actual  deposition  in  the  channel  for 
the  two  channel  configurations. 


The  bed  shear  stress  is  estimated  using  the  quadratic  relationship 
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T  h=\Cf9V2  (2) 

where: 

Cf  =  a  coefficient  of  friction 
p  =  density 
V  =  velocity 

The  value  of  C  f  depends  upon  the  nature  of  the  velocity  being  used.  If  the 

velocity  is  the  near-bed  velocity  the  coefficient  can  be  large.  In  this  case, 
the  velocity  is  a  cross-section  average  value  and  the  coefficient  is  relatively 
small.  For  established  open  channel  flow,  the  coefficient  is  about  0.005 
(Maynord  2005b).  However,  in  this  case,  the  boundary  layer  is  quite  thin 
since  the  flow  profile  does  not  have  time  to  reach  a  quasi-steady  condition. 
In  this  case  Maynord  (2005b)  recommends  using  Cf  =0.01.  Using  this 

coefficient  one  now  estimates  the  shear  stress  for  the  largest  vessel  in  the 
two  channel  configurations  (see  Table  4). 


Table  4.  Erosion  rate  for  a  vessel  passage. 


Channel  Condition 

Return  Current  (m/s) 

(Vr) 

Shear  Stress  (Pa) 

(TJ 

Erosion  per 
passage  (m)  ( ds ) 

Base 

1.3 

8.5 

8.8xl0-6 

Plan 

1.1 

6.1 

4.6x10-6 

The  coefficients  needed  to  calculate  the  impact  upon  sedimentation  are 
derived  from  the  field  sediment  data  collection.  Figure  17  shows  the  ero¬ 
sion  rate  compared  to  the  shear  stress  at  push  core  sta  1  in  the  bay  west  of 
the  Houston  Ship  Channel.  From  inspection,  the  critical  shear  stress  for 
erosion,  xce  (see  Partheniades  1962),  is  about  1.0  Pa. 

Kg 

The  erosion  rate,  Ke ,  is  calculated  to  be  6.8x10  — - - - - from  this 

nrx  s  x  Pa 

figure  as  the  slope  of  the  line  between  0.97  Pa  and  1.14  Pa.  Assume  that 
the  dry  sediment  density,  ps ,  is  500  kg/m3.  The  depth  of  erosion  per  ves¬ 
sel  passage,  ds ,  is  calculated  using  the  following  equation. 

'p 

ds  =77  JX(u>-Ue)* 

0 


(3) 
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where 


xb  =  bed  stress 
tp  =  time  of  passage 


This  time  of  passage,  tp ,  is  the  length  of  time  the  bed  is  subject  to  the 

stress.  In  this  case  it  is  assumed  to  be  the  length  of  the  vessel  divided  by 
the  speed  of  the  vessel.  All  the  terms  are  assumed  to  be  a  constant.  The 
results  are  shown  in  the  last  column  of  Table  4.  From  the  2004  records  of 
the  Port  of  Houston  (Dept,  of  the  Army  2006)  there  were  4280  passages  of 
vessels  with  a  draft  of  30  ft  or  more  in  this  year.  Using  4280  passages  in  a 
year  the  total  erosion  for  base  conditions  would  be  about  0.04  meters  per 
year  and  for  the  plan  condition  about  0.02  meters  per  year.  The  dredging 
records  suggest  that  channel  deposition  in  the  upper  bay  is  on  the  order  of 
0.4  meters  per  year.  This  would  suggest  that  the  increase  in  apparent 
shoaling  due  to  the  reduced  stress  of  the  return  currents  associated  with 
the  vessel  passage  is  on  the  order  of  5%.  This  does  not  include  the  fact  that 
larger  vessels  will  utilize  the  larger  channel  thus  making  this  effect  even 
less.  The  analysis  does  not  include  erosion  caused  by  the  ship  propeller.  It 
has  been  noted  in  the  field  data  section  that  the  center  of  the  channel  has 
maintained  its  “as  built”  condition.  The  channel  edges  however  have 
shoaled  considerably.  This  is  most  likely  a  result  of  the  vessel  propeller 
with  some  influence  from  the  vessel  return  currents.  Our  estimate  of  the 
erosion  caused  by  ship  passage  alone  only  accounts  for  about  0.02  meters 
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of  increased  deposition  annually.  Therefore,  the  difference  in  vessel  ero¬ 
sion  in  the  channel  with  the  two  different  channel  sizes  doesn't  appear  to 
be  large  enough  to  explain  the  deposition  rate  shown  in  the  enlarged 
channel. 
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5  Numerical  Model  Simulation  and  Analysis 

This  chapter  details  the  model  study  analysis  of  the  hydrodynamics  and  a 
sediment  tracer  in  Galveston  Bay,  TX.  TABS-MDS  is  the  finite  element 
numerical  model  used  for  this  investigation,  which  models  three- 
dimensional  (3-D)  hydrodynamics  and  salinity  as  well  as  suspended  sedi¬ 
ment  transport  and  bed  changes.  This  model  has  been  widely  used  by 
ERDC  to  model  3-D  hydrodynamics  and  salinity  at  numerous  locations, 
including  Galveston  Bay  (Berger  et  al.  1995a, b).  Numerical  model  simula¬ 
tions  were  performed  for  both  the  original  40-  x  400-ft  (12.2-  x  122-m) 
channel  and  the  45-  x  530-ft  (13.7-  x  162-m)  channel.  The  enlarged  chan¬ 
nel  dimensions  are  nominally  45-  x  530-ft  (13.7-  x  162-m);  however,  actual 
surveys  were  used  and  the  depth  for  the  45-ft-  (13.7-m-)  deep  channel  is 
typically  about  49  ft  (14.9  m)  mlw  at  the  channel  centerline.  Initial  runs 
were  made  simulating  hydrodynamic  effects  only  in  order  to  determine 
residual  shear  stresses  on  the  bed.  These  stresses  give  an  indication  of  any 
changes  in  the  bed  due  to  velocity  variations  among  the  two  channel  con¬ 
figurations.  Sediment  conditions  were  then  added  to  the  simulations. 
These  runs  do  not  represent  a  full  sediment  simulation  of  the  model  do¬ 
main,  but  rather  give  an  estimate  of  where  suspended  sediment  travels 
and  settles  given  the  grain  characteristics,  model  hydrodynamics,  and 
wind  induced  wave  effects. 

Experimental  Conditions 

The  model  conditions  use  the  present  low-flow  and  medium-flow  hydro- 
logic  eras  from  the  previous  investigations.  The  river  inflows  for  the  Trin¬ 
ity  River  and  San  Jacinto  River  were  included  in  the  model  conditions. 

The  tidal  boundary  data  were  from  NOS/NOAA  sta  877-1510  at  Pleasure 
Pier  for  the  calendar  year  1984,  which  was  a  typical  temperature  year.  The 
model  boundary,  located  26  miles  (42  km)  offshore,  was  applied  by  shift¬ 
ing  these  data  in  time  by  1.31  hr  to  match  the  phase  at  Pleasure  Pier.  Peri¬ 
ods  of  less  than  3  hr  were  removed  by  filtering  in  order  to  prevent  noise  in 
the  results.  Wind  data  from  1984  were  obtained  from  the  National 
Weather  Service  for  Houston  Intercontinental  Airport.  These  data  were 
transposed  to  the  bay  system  after  a  correlation  with  field  data  collected 
for  the  prior  studies.  Gulf  boundary  salinities  are  based  on  15-year 
monthly  averages  for  the  Gulf  and  the  initial  salinity  field  was  taken  from 
the  final  salinity  field  produced  in  the  verification  runs.  For  a  detailed  de¬ 
scription  of  the  experimental  conditions,  see  Berger  et  al.  1995b. 
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Analysis  of  Model  Currents 

This  section  describes  the  residual  shear  stresses  in  Galveston  Bay  in 
and  around  the  Houston  Ship  Channel  for  the  conditions  of  the  40-ft- 
(12.2-m-)  deep  by  400-ft-  (122-m-)  wide  channel  and  for  the  new  channel 
configuration  of  45-ft-  (13.7-m-)  deep  by  530-ft-  (162-m-)  wide.  It  also  es¬ 
timates  a  shoaling  volume  for  the  two  channel  configurations.  The  40-  x 
400-ft  (12.2-  x  122-m)  channel  is  termed  the  "base"  conditions  and  the  45- 
x  530-ft  (13.7-  x  162-m),  the  “plan”  conditions.  The  40-  x  400-ft  (12.2-  x 
122-m)  channel  existed  in  Galveston  Bay  from  the  early  1960's  into  the  late 
1990's.  At  that  time  the  channel  was  enlarged  to  the  plan  conditions.  The 
enlarged  channel  has  shown  a  significant  increase  in  required  dredging. 
Planners  need  to  know  if  this  new  increased  dredging  requirement  is  a 
permanent  feature  or  is  the  shoaling  going  to  return  to  the  earlier  channel 
dredging  rates.  Typically  the  dredging  rate  is  elevated  while  new  construc¬ 
tion  is  taking  place.  This  results  from  leakage  from  the  dredging  operation 
and  placement  as  well  as  the  channel  shoulder  instability.  As  a  part  of  that 
evaluation,  the  currents  from  the  3-D  model  of  the  bay  are  analyzed  to  de¬ 
termine  if  the  currents  have  changed  enough  to  produce  the  increase  in 
sedimentation. 

An  indicator  of  the  average  shear  stress  (the  velocity  squared)  is  plotted 
for  a  few  representative  months,  made  for  the  near-bed  conditions  only.  A 
location  map  for  use  in  this  chapter  is  shown  in  Figure  18. 

Monthly  Shear  Residuals 

These  residual  shear  stresses  are  given  by  month.  The  shear  stress  values 
indicate  the  direction  that  material  in  periodic  contact  with  the  bed  would 
tend  to  move.  They  also  reflect  the  direction  that  suspended  sediment 
moves  and  are  offered  as  a  means  to  understand  the  major  hydrodynamic 
features  of  the  bay  sediment  movement  before  launching  a  full  sediment 
model. 

These  shear  residuals  are  actually  the  mean  velocity  squared  for  each 
month.  This  is  calculated  as  follows: 

T 


j  (u(t)i  +  v(/)j)|V(/)|t// 


where: 
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V  =  w(t)  +  v(t)j 

T  =  time  period  over  which  the  residual  is  being  evaluated. 


Figure  18.  Location  map  of  Galveston  Bay,  showing  depth  contours  and  velocity 

point  data  locations. 


The  runs  used  in  the  analysis  were  yearlong  simulations.  However,  the 
first  few  months  are  used  to  let  the  system  adjust  from  initial  conditions. 
The  salinity  is  slower  to  acclimate  than  are  the  water-surface  elevations 
and  currents.  Allowing  the  first  3  months  for  this  spin-up  period,  and 
since  the  months  of  November  and  December  are  much  like  October,  only 
April  through  October  are  included.  Figures  19  through  32  contain  the 
midbay  residual  shear  results  and  Figures  33  through  46  show  those  for 
the  upper  bay. 

First  consider  the  midbay  shear  residual  for  the  40-  x  400-ft  (12.2-  x 
122-m)  channel  in  Figures  19-25.  April  shows  a  strong  flood  residual 
shear.  This  weakens  in  May,  but  is  still  flood-oriented.  June  and  July 
have  strong  ebb  residual  shear  down  to  about  Redfish  Reef,  where  the  re¬ 
sidual  shear  is  very  small.  One  might  expect  sediment  accumulation  in 
this  region,  at  this  time.  August  through  October  show  a  return  to  flood- 
directed  shear. 
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The  midbay  shear  residual  for  the  45-  x  530-ft  (13.7-  x  162-m)  channel  is 
shown  in  Figures  26-32.  All  months  show  flood  oriented  shear  residuals. 
However,  they  are  weak  in  June  and  are  strong  in  April,  September,  and 
October.  There  is  also  some  noticeable  shear  shown  in  the  shallows  during 
September  and  October  around  the  midbay  marsh  site. 

In  the  upper  bay,  the  shear  residual  for  the  40-  x  400-ft  (12.2-  x  122-m) 
channel  is  shown  in  Figures  33-39.  The  April  pattern  is  for  a  strong  flood 
residual  in  the  channel,  except  near  the  Bayport  Channel  flare  where  the 
residual  is  essentially  zero.  The  shallows  near  the  channel  tend  to  have  an 
ebb  residual  shear.  The  May  pattern  is  a  weak  version  of  April.  June  and 
July  have  ebb  residual  shear  in  the  channel.  There  is  also  a  flood  shear  re¬ 
sidual  on  the  western  edge  of  the  bay  in  the  shallows.  In  the  subsequent 
months  the  pattern  returns  to  that  of  April,  though  weakly  in  August. 

The  upper  bay  shear  residual  for  the  45-  x  530-ft  (13.7-  x  162-m)  channel  is 
shown  in  Figures  40-46.  The  April  pattern  shows  a  strong  flood  residual 
in  the  channel,  a  dead  area  near  the  Bayport  Channel  flare,  and  an  ebb  re¬ 
sidual  in  the  shallows  near  the  channel.  This  pattern  is  similar  for  the 
other  months,  but  is  weak  during  June  and  strengthens  in  later  months. 

Figure  47  shows  the  total  freshwater  inflow  to  Galveston  Bay  in  these 
model  tests.  Figure  48  shows  the  salinity  used  at  the  Gulf  boundary  dur¬ 
ing  these  runs.  A  more  in-depth  explanation  for  these  particular  boundary 
conditions  is  found  in  Berger  et  al.  (1995b).  This  is  a  simulated  low  fresh¬ 
water  flow  year.  Therefore,  the  saltwater  is  not  being  pushed  out  of  the 
bay.  In  fact,  the  salinity  in  the  upper  Galveston  Bay  is  around  15  ppt  dur¬ 
ing  the  June-July  period  for  the  40-  x  400-ft  (12.2-  x  122-m)  channel. 

What  appears  to  actually  be  the  cause  for  the  difference  of  the  June  and 
July  shear  directions  for  both  the  channel  and  the  shallows  is  the  drop  in 
Gulf  salinity  to  a  minimum  at  mid-May.  This  drop  in  salinity  results  in  an 
evacuation  of  saline  water  from  the  bay.  The  salt  water  is  retreating  from 
the  deeper  areas,  such  as  the  channel.  When  this  salt  water  is  exiting  the 
bay,  it  is  being  replaced  with  fresher  water  moving  upland  in  the  shallows. 
This  is  reversed  from  the  normal  pattern  of  ebb  directed  shears  in  the  shal¬ 
lows  and  flood  directed  net  shears  in  the  channel. 

Individual  Shear  Locations 

The  results  at  eight  individual  locations  throughout  the  bay  and  channel 
are  shown  in  Figures  49-72.  Locations  1  through  6  are  distributed  in  the 
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upper  bay.  Locations  lC  and  2C  are  in  the  channel.  Each  figure  shows  the 
percentage  of  time  for  which  the  velocity  squared  is  exceeded.  While  the 
solutions  were  calculated  over  a  year,  the  months  of  April,  June,  and  Octo¬ 
ber  were  used  as  indicators  and  they  are  shown  in  these  figures.  For  the 
months  of  April  and  June  the  currents  were  generally  reduced  as  a  result 
of  the  channel  enlargement.  This  is  particularly  true  in  the  channel.  It  is 
also  noted  that  at  the  stations  farther  from  the  channel  (Points  2,  4,  and  6) 
the  reduction  in  currents  is  more  than  the  near  channel  stations  (Points  1, 
3,  and  5).  It  is  important  to  note  that  for  April  and  June,  the  near-bed 
channel  currents  have  been  reduced.  This  could  lead  to  increased  deposi¬ 
tion.  For  all  points,  the  October  currents  did  not  vary  greatly  between  the 
two  channel  configurations. 

Generally,  the  Houston  Ship  Channel  within  Galveston  Bay  is  flood- 
directed  shear.  However,  for  the  base  conditions  in  the  months  of  June 
and  July,  the  shear  is  ebb-directed  in  the  northern  portion  of  the  bay  chan¬ 
nel.  This  period  of  ebb-dominated  shear  is  likely  due  to  the  drop  in  salin¬ 
ity  in  the  model  ocean.  The  dense  water  is  then  tending  to  evacuate  the 
system.  The  period  of  low  salinity  at  the  Gulf  of  Mexico  is  due  to  the 
freshwater  currents  from  the  Atchafalaya  and  Mississippi  Rivers.  For  the 
other  months  and  for  the  plan  channel  the  bed  shear  is  flood-directed.  In 
June  and  July  the  plan  condition  shears  are  flood-directed,  but  much  less 
strongly. 

The  area  around  the  junction  of  the  Bayport  Channel  and  the  Houston 
Ship  Channel  shows  a  disrupted  shear  pattern.  The  shear  is  low  in  the 
flare  for  the  Bayport  Channel  but  also  in  the  Houston  Ship  Channel  in  that 
location.  This  does  not  appear  when  the  shear  is  ebb-directed,  but  it  does 
whenever  the  shear  is  flood-directed.  There  is  a  bend  in  the  Houston  Ship 
Channel  just  south  of  this  junction  that  may  be  the  cause. 

The  shear  at  several  individual  locations  is  presented.  Of  particular  inter¬ 
est  are  two  channel  stations  that  show  that  the  enlarged  channel  results  in 
reduced  shear.  This  reduction  in  shear  should  result  in  an  increase  in 
deposition. 
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Figure  19.  April,  midbay,  bottom  depth,  40X400  channel,  residual  shear. 


Figure  20.  May,  midbay,  bottom  depth,  40X400  channel,  residual  shear. 
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Figure  21.  June,  midbay,  bottom  depth,  40X400  channel,  residual  shear. 


Figure  22.  July,  midbay,  bottom  depth,  40X400  channel,  residual  shear. 


ERDC/CHL  TR-06-8 


38 


Figure  23.  August,  midbay,  bottom  depth,  40X400  channel,  residual  shear. 


Figure  24.  September,  midbay,  bottom  depth,  40X400  channel,  residual  shear. 
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Figure  25.  October,  midbay,  bottom  depth,  40X400  channel,  residual  shear. 


Figure  26.  April,  midbay,  bottom  depth,  45X530  channel,  residual  shear. 
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Figure  27.  May,  midbay,  bottom  depth,  45X530  channel,  residual  shear. 


Figure  28.  June,  midbay,  bottom  depth,  45X530  channel,  residual  shear. 
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Figure  29.  July,  midbay,  bottom  depth,  45X530  channel,  residual  shear. 


Figure  30.  August,  midbay,  bottom  depth,  45X530  channel,  residual  shear. 
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Figure  31.  September,  midbay,  bottom  depth,  45X530  channel,  residual  shear. 


Figure  32.  October,  midbay,  bottom  depth,  45X530  channel,  residual  shear. 
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Figure  33.  April,  upper  bay,  bottom  depth,  40X400  channel,  residual  shear. 


Figure  34.  May,  upper  bay,  bottom  depth,  40X400  channel,  residual  shear. 
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Figure  35.  June,  upper  bay,  bottom  depth,  40X400  channel,  residual  shear. 


Figure  36.  July,  upper  bay,  bottom  depth,  40X400  channel,  residual  shear. 
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Figure  37.  August,  upper  bay,  bottom  depth,  40X400  channel,  residual  shear. 


Figure  38.  September,  upper  bay,  bottom  depth,  40X400  channel,  residual  shear. 
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Figure  39.  October,  upper  bay,  bottom  depth,  40X400  channel,  residual  shear. 


Figure  40.  April,  upper  bay,  bottom  depth,  45X530  channel,  residual  shear. 
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Figure  41.  May,  upper  bay,  bottom  depth,  45X530  channel,  residual  shear. 


Figure  42.  June,  upper  bay,  bottom  depth,  45X530  channel,  residual  shear. 
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Figure  43.  July,  upper  bay,  bottom  depth,  45X530  channel,  residual  shear. 


Figure  44.  August,  upper  bay,  bottom  depth,  45X530  channel,  residual  shear. 
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Figure  45.  September,  upper  bay,  bottom  depth,  45X530  channel,  residual  shear. 


Figure  46.  October,  upper  bay,  bottom  depth,  45X530  channel,  residual  shear. 


ERDC/CHL  TR-06-8 


50 


Figure  47.  Flow  rate  in  combined  Trinity  and  San  Jacinto  rivers  for  low  flow, 

1990  hydrology  year. 


Figure  48.  Salinity  used  for  Gulf  boundary  in  model  simulations. 
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Figure  49.  April  velocity  squared  exceedance  for  Point  1. 


Figure  50.  April  velocity  squared  exceedance  for  Point  2. 
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Figure  51.  April  velocity  squared  exceedance  for  Point  3. 


Figure  52.  April  velocity  squared  exceedance  for  Point  4. 
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Figure  53.  April  velocity  squared  exceedance  for  Point  5. 


Figure  54.  April  velocity  squared  exceedance  for  Point  6. 
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Figure  56.  April  velocity  squared  for  Point  2C. 
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Figure  57.  June  velocity  squared  exceedance  for  Point  1. 


Figure  58.  June  velocity  squared  exceedance  for  Point  2. 
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Figure  59.  June  velocity  squared  exceedance  for  Point  3. 


Figure  60.  June  velocity  squared  exceedance  for  Point  4. 
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Figure  61.  June  velocity  squared  exceedance  for  Point  5. 


Figure  62.  June  velocity  squared  exceedance  for  Point  6. 
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Figure  63.  June  velocity  squared  exceedance  for  Point  1C. 


Figure  64.  June  velocity  squared  exceedance  for  Point  2C. 
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Figure  66.  October  velocity  exceedance  for  Point  2. 
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Figure  67.  October  velocity  squared  exceedance  for  Point  3. 


Figure  68.  October  velocity  squared  exceedance  for  Point  4. 
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Figure  69.  October  velocity  squared  exceedance  for  Point  5. 


Figure  70.  October  velocity  squared  exceedance  for  Point  6. 


ERDC/CHL  TR-06-8 


62 


Figure  71.  October  velocity  squared  exceedance  for  Point  1C. 


Figure  72.  October  velocity  squared  exceedance  for  Point  2C. 
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Sediment  Estimates  based  upon  Hydrodynamic  Model 

As  a  preliminary  effort  the  shears  obtained  from  the  hydrodynamic  model 
are  used  to  estimate  the  shoaling  rate  in  the  channel.  This  part  of  the  re¬ 
port  was  done  before  the  sediment  model  was  in  operation  and  was  used 
to  provide  a  quick  estimate.  This  is  followed  up  in  the  next  chapter  with 
the  sediment  model,  at  which  time  some  of  the  parameters  used  in  this 
section  may  be  modified.  A  rough  estimate  of  the  shoaling  rate  during 
April  can  be  made  using  values  from  the  field  data  and  other  bounding  es¬ 
timates.  This  estimate  will  allow  us  to  more  precisely  define  the  critical 
shear  stress  for  deposition.  This,  in  turn,  will  be  used  by  the  model  to  cal¬ 
culate  the  deposition  rate  change  caused  by  the  channel  enlargement. 

April  in  the  field  is  typically  a  high  shoaling  period.  An  estimate  for  the 
coefficient  of  friction  is  Cf  =  0.091  in  the  equation  x  =  ■|pC/V'2 .  This  es¬ 
timate  comes  from  the  Manning  coefficient  used  in  the  model  run  of  0.033 
for  the  channel.  This  friction  is  based  on  the  bottom  velocity  and  so  the 
Manning  coefficient  is  higher  than  when  using  the  mean  velocity.  The  ve¬ 
locity-squared  is  in  feet  per  second  and  shear  stress  is  calculated  in  Pascals 
(Pa).  The  rate  of  deposition  can  be  estimated  as 

d  =  VsPCB  (Krone  1962)  (5) 


where: 


V. 


P 


CB 


=  settling  velocity 

=  probability.  That  reduces  the  settling  velocity  to  an 

effective  settling  velocity;  if  x  >  xcd,P  =  0 
=  bed  concentration. 


The  estimates  of  the  shear  exceedance  for  locations  lC  and  2C  are  shown 
in  Figures  55  and  56,  respectively.  The  following  estimated  parameters 
were  found  using  representative  values  from  the  collected  field  data  of  Au¬ 
gust  2004. 

CB  =  200  mg/L  (0.2  kg/m3) 

Vs  =  0.5  mm/s  (0.0005  m/s) 


Other  estimated  parameters  are  as  follows: 
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Ss  =  2.65  is  the  grain  specific  gravity  commonly  used  value 

Dry  sediment  density  =  500  kg/m3 

Critical  shear  stress  for  deposition  rcd  -0.5  Pa. 

These  final  two  parameters  were  suggested  by  Teeter  (2005).  The  dry 
sediment  density  is  for  fresh  deposits.  Over  time,  the  deposit  will  consoli¬ 
date  and  the  density  will  increase.  The  field  data  found  that  the  critical 
shear  stress  for  erosion  in  the  bay  was  about  1.0  to  1.2  Pa.  This  would  be 
the  upper  limit  for  xcd .  Also  samples  in  this  area  of  the  channel  contained 

fine  sand  and  medium  to  coarse  silts  (along  with  some  clay).  This  suggests 
a  value  of  approximately  0.5  Pa.  Values  of  0.25,  0.5,  and  0.75  were  used  in 
the  calculations  simply  to  test  sensitivity. 

Tables  5  and  6  contain  estimates  of  deposition  for  April  for  Point  lC  and 
2C,  respectively.  The  first  major  heading  is  "Probability"  which  is  the  inte¬ 
gral  of  the  instantaneous  probability  over  the  month.  The  second  major 
heading  is  the  deposition  rate  in  meters/month.  And  the  final  major  head¬ 
ing  is  the  volume  deposition  per  meter  of  channel  length.  This  column  is 
found  by  multiplying  the  deposition  by  the  channel  width.  This  width  is 
400  ft  (122  m)  for  the  base  conditions. 

The  different  critical  shear  stresses  give  typical  depth  deposit  rates  for  the 
base  conditions  between  0.10  and  0.38  m/month  for  Point  lC  and  between 
0.06  and  0.36  for  Point  2C.  The  lower  limit  of  deposition  is  for  xcd  = 

0.25  Pa.  The  deposition  amounts  for  this  critical  shear  stress  for  deposi¬ 
tion  are  too  small  and  the  relative  size  of  deposition  at  the  two  points  is 
strongly  weighted  toward  Point  lC.  The  actual  dredging  records  show  that 
the  two  points  have  similar  shoaling  rates.  For  the  critical  shear  stress  for 
deposition  of  0.5  and  0.75  Pa,  the  distribution  is  more  even  and  the  rate  is 
more  like  what  has  been  observed  in  the  field.  Therefore,  the  results  for 
0.5  Pa  are  probably  the  most  reliable  and  will  be  used  in  subsequent  calcu¬ 
lations. 


Table  5.  Estimates  of  shoaling  rate  around  Point  1C. 


Point  1C 

April 

Ted ,  Pa 

Probability 

Deposition  Rate, 
m/month 

Volume  (m3)  per  meter 
length 

Base 

Base 

Base 

0.25 

0.18 

0.10 

12 

0.50 

0.59 

0.31 

38 
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0.75 

0.73 

0.38 

46 

Table  6.  Estimates  of  shoaling  rate  around  Point  2C. 


Point  2C 

April 

Ted  Pa 

Probability 

Deposition  Rate, 
m/month 

Volume  (m3)  per  meter 
length 

Base 

Base 

Base 

0.25 

0.12 

0.06 

8 

0.50 

0.54 

0.28 

34 

0.75 

0.69 

0.36 
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Sediment  Model  Results 

This  chapter  utilizes  a  sediment  model  in  which  a  tracer  sediment,  that  can 
settle  and  erode  based  upon  the  shear  stress,  is  introduced  to  the  system  to 
estimate  the  impact  of  the  channel  enlargement.  In  this  chapter  the  sedi¬ 
ment  model  is  used  in  an  unvalidated  state.  The  next  study  (Phase  2)  will 
validate  the  model  and  use  it  to  evaluate  alternatives  to  reduce  the  re¬ 
quired  channel  maintenance.  This  chapter  is  a  follow  on  of  the  previous 
chapters.  These  chapters  were  provided  to  the  sponsor  as  they  were  com¬ 
pleted  to  provide  estimates  as  early  as  possible.  The  adjustment  period  in 
this  model  began  with  the  parameters  developed  in  the  prior  chapters  but 
in  order  to  behave  more  like  the  system  these  were  adjusted  as  needed. 
This  model  study  made  no  attempt  to  introduce  the  sediment  via  the  rivers 
or  attempt  to  match  quantities  of  shoaling.  Here  the  model  used  a  tracer 
sediment  introduced  along  the  channel.  The  ultimate  deposition  distribu¬ 
tion  is  compared  with  the  field  distribution.  This  sediment  tracer  is  then 
introduced  over  the  same  area  with  the  enlarged  channel  condition.  The 
change  in  deposition  distribution  and  amount  is  then  compared  with  the 
base  (40X400  ft)  channel  condition. 

Hydrodynamic  and  salinity  coupled  runs  were  made  using  the  "Present" 
low-flow  and  medium-flow  hydrologic  eras  from  a  prior  study  as  described 
earlier  (see  Berger  et  al.  1995b).  The  effects  of  wind  waves  on  the  sediment 
movement  were  included  in  the  TABS-MDS  simulations.  The  wind  was 
applied  using  the  same  speeds  and  directions  as  in  the  hydrodynamic 
simulations.  These  runs  were  made  from  January  through  June.  A  single 
grain  size  sediment  was  “rained”  into  the  model  along  the  channel  region 
for  both  configurations.  An  identical  sediment  load  was  added  to  both 
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domains  over  the  same  areas.  The  sediment,  with  a  concentration  of  l 
kg/m3,  was  introduced  for  3  days  at  a  rate  of  0.01  m/s  over  the  area  shown 
in  Figure  73.  The  sediment  was  modeled  for  April,  May,  and  June  sepa¬ 
rately  with  grain  and  bed  characteristics  determined  from  the  field  data 
collection.  The  bed  was  defined  such  that  the  only  sediment  available  for 
erosion  was  that  which  entered  during  the  run.  Each  run  was  started  10 
hours  prior  to  midnight  on  the  first  of  the  month  in  order  to  allow  the 
model  to  spin  up  before  introducing  the  sediment.  The  characteristics  of 
the  sediment  used  in  the  model  are  provided  in  Table  7. 


Table  7.  Sediment  characteristics 


Sediment  Characteristics 

Density 

1.565  g/cm3 

Erosion  Rate 

0.0000346  kg/m2*sec 

Critical  Shear  for  Erosion 

1.336  Pa 

Critical  Shear  for  Deposition  (Tcd)* 

0.5  Pa 

Settling  Velocity  ( Vs  )* 

0.5  mm/sec 

*  denotes  variables  modified  for  sensitivity  analysis 

The  ratios  of  the  volumetric  change  in  shoaling  along  the  channel  reach 
between  the  two  configurations  for  both  the  low-flow  and  medium-flow 
are  given  in  Table  8.  These  are  given  separately  for  each  of  the  three 
months  simulated. 


Table  8.  Shoaling  volume  ratios  for  low  and  medium  flows. 


Shoaling  Volume  Ratio  (Plan/Base) 

Flow  Condition 

April 

May 

June 

Low 

1.28 

1.19 

1.17 

Medium 

1.34 

1.18 

1.07 

The  increase  in  volume  of  deposition  in  the  channel  ranged  from  7%  in 
June  for  the  medium  flow  event  to  34%  in  April  for  the  same  flow  condi¬ 
tion.  Given  the  limitations  of  these  sediment  simulations,  the  predicted 
increase  in  shoaling  within  the  enlarged  channel  could  conservatively  be 
predicted  to  range  from  20-30  percent. 
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Figure  73.  40-  X  400-ft  channel  configuration  and  location  map,  blue  dots  are  3,000-m 
increments  from  0  at  Tabb’s  Bay  to  42,000  m  at  Bolivar  Roads. 


Figures  74-76  show  the  change  in  distribution  between  the  base  condition 
and  the  plan  (enlarged  channel)  condition  for  the  low-flow  and  medium- 
flow  hydrology.  These  figures  only  consider  the  Houston  Ship  Channel 
and  do  not  account  for  any  shoaling  outside  of  the  channel.  The  distribu¬ 
tion  plots  are  given  in  3,281-ft  (1,000-m)  increments  that  can  be  refer¬ 
enced  using  the  9,842.5-ft  (3,000-m)  markers  (blue  dots)  in  Figure  73. 
Generally,  the  distribution  has  shifted  to  show  more  shoaling  toward  the 
upland  regions.  In  April  there  is  also  a  shift  toward  the  Red  Fish  Reef 
channel  region. 
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45X530  ft  Percentage  Increase  in  Sediment  Volume  from  40X400  ft 

May 


Distance  from  Tabbs  Bay  (1000  meters) 


□  Low  Flow  ■  Medium  Flow 


Figure  75.  Percentage  increase  in  shoaling  for  enlarged  channel  for  May. 
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45X530  ft  Percentage  Increase  in  Sediment  Volume  from  40X400  ft 

June 


Distance  from  Tabbs  Bay  (1000  meters) 


□  Low  Flow  ■  Medium  Flow 


Figure  76.  Percentage  increase  in  shoaling  for  enlarged  channel  for  June. 

Sensitivity  analyses  were  performed  to  determine  which  variables  affected 
the  shoaling  results  given  by  the  sediment  simulations.  The  settling  veloc¬ 
ity  ( 1/  )  and  critical  shear  stress  for  deposition  ( xcd )  were  varied  as  a 

means  to  influence  the  time  particles  remained  in  the  water  column.  The 
domain  is  generally  a  depositional  one  and  so  the  critical  shear  stress  for 
erosion  has  little  effect.  Lowering  the  critical  shear  stress  for  deposition 
will  allow  the  sediment  to  deposit  at  lower  velocities.  Lowering  the  set¬ 
tling  velocity  will  keep  the  sediment  in  suspension  longer  and  allow  it  to 
move  further  throughout  the  system.  So  by  lowering  both  of  these  values 
independently,  the  possibility  of  high  and  low  deposition  throughout  the 
length  of  the  channel  will  be  included  in  the  sensitivity  analyses.  Both  of 
these  values  were  varied  from  a  magnitude  of  0.5  to  0.1  in  order  to  deter¬ 
mine  the  impact  of  extreme  variation  in  these  parameters.  These  varia¬ 
tions  did  show  a  change  in  the  volume  of  the  deposited  sediment  within 
the  channel  among  the  configurations;  however,  this  change  was  not  large 
and  the  distribution  of  the  sediment  did  not  vary.  The  wind  wave  effects 
were  included  and  removed  as  well.  In  this  case,  the  wind  waves  would 
only  be  able  to  reintroduce  sediment  that  had  fallen  in  the  shallows  during 
the  simulation  since  the  initial  state  of  the  model  was  to  have  no  sediment 
on  the  bed.  Tables  9-16  give  the  ratios  by  which  the  shoaling  volume 
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changed  for  the  base  and  sensitivity  runs  for  both  flow  conditions  within 
the  channel  from  Tabb’s  Bay  to  Bolivar  Roads,  within  the  channel  along 
Atkinson  Island,  in  the  Bayport  Flare  and  Channel,  and  in  Barbour’s  Cut. 
The  Bayport  and  Barbour's  Cut  channels  were  not  enlarged  in  this  study. 
The  effects  of  enlargement  are  only  from  the  Houston  Ship  Channel.  The 
original  conclusion  of  a  20  to  30%  increase  seems  reasonable  even  consid¬ 
ering  the  impact  of  this  sensitivity. 


Table  9.  Shoaling  volume  ratios  for  channel  and  low  flow. 


Low  Flow  Hydrology 

Channel  Only 

Change  Factor 

April 

May 

June 

40  X  400  Base 

1.28 

1.19 

1.17 

vs  Tcd  =0.1  Pa 

1.34 

1.02 

0.99 

45X530  Vs  =0.1  mm/s 

1.22 

1.21 

1.23 

Table  10.  Shoaling  volume  ratios  for  channel  and  medium  flow. 


Medium  Flow  Hydrology 

Channel  Only  1 

Change  Factor 

April 

May 

June 

40X400  Base 

1.34 

1.18 

1.07 

vs  Tcd  =0.1  Pa 

1.54 

0.99 

0.91 

45X530  Vs  =0.1  mm/s 

1.31 

1.21 

1.19 

Table  11.  Shoaling  volume  ratios  for  Atkinson  Island  and  low  flow. 


Low  Flow  Hydrology 

Atkinson  Island 

Change  Factor 

April 

May 

June 

40X400  Base 

1.26 

1.20 

1.20 

vs  Tcd  =0.1  Pa 

1.24 

1.10 

1.04 

45X530  Vs  =0.1  mm/s 

1.12 

1.18 

1.27 

Table  12.  Shoaling  volume  ratios  for  Atkinson  Island  and  medium-flow. 


Medium  Flow  Hydrology 

Atkinson  Island 

Change  Factor 

April 

May 

June 

40X400  Base 

1.36 

1.22 

1.19 

vs  7 cd  =  0.1  Pa 

1.56 

1.11 

0.98 

45X530  Vs  =0.1  mm/s 

1.24 

1.26 

1.16 

Table  13.  Shoaling  volume  ratios  for  Bayport  Channel  and  Flare  and  low  flow. 


Low  Flow  Hydrology 

Bayport  Flare 

Change  Factor 

April 

May 

June 

40X400  Base 

1.08 

0.99 

1.07 

vs  Tcd  =0.1  Pa 

1.09 

0.94 

1.07 

45X530  Vs  =0.1  mm/s 

0.93 

1.13 

1.21 
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Table  14.  Shoaling  volume  ratios  for  Bayport  Channel  and  Flare  and  medium  flow. 


Medium  Flow  Hydrology 

Bayport  Flare  1 

Change  Factor 

April 

May 

June 

40X400  Base 

1.01 

1.05 

1.00 

vs  Tcd  =0.1  Pa 

0.96 

1.02 

0.92 

45X530  Vs  =0.1  mm/s 

1.01 

1.12 

1.04 

Table  15.  Shoaling  volume  ratios  for  Barbour’s  Cut  and  low  flow. 


Low  Flow  Hydrology 

Barbour's  Cut 

Change  Factor 

April 

May 

June 

40X400  Base 

0.98 

0.99 

1.04 

vs  Tcd  =  0.1  Pa 

1.00 

1.01 

1.05 

45X530  Vs  =0.1  mm/s 

0.95 

0.97 

1.08 

Table  16.  Shoaling  volume  ratios  for  Barbour’s  Cut  and  medium  flow. 


Medium  Flow  Hydrology 

Barbour's  Cut 

Change  Factor 

April 

May 

June 

40X400  Base 

1.19 

1.02 

1.18 

vs  Tcd  =0.1  Pa 

1.21 

1.03 

1.21 

45X530  Vs  =0.1  mm/s 

1.27 

1.07 

1.44 
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6  Summary  and  Conclusions 

This  report  analyzes  the  likely  causes  of  an  increase  in  channel  sedimenta¬ 
tion  due  to  the  enlargement.  It  has  been  noted  that  a  large  increase  has 
occurred  since  the  enlarged  channel  has  been  in  place  and  the  question  is 
whether  or  not  this  is  a  permanent  increase  or  will  the  sedimentation  rate 
drop  back  to  more  typical  levels.  This  investigation  takes  place  and  is  re¬ 
ported  in  stages.  As  each  stage  was  completed  the  results  were  reported  to 
the  sponsor.  This  report  is  a  compilation  of  each  of  these  interim  reports. 

Dredging  Analysis 

The  historic  dredging  records  indicate  that  the  shoaling  rate  was  higher  in 
the  upper  half  of  the  bay  channel  than  in  the  lower  bay  channel  for  the  40- 
x  400-ft  (12.2-  x  122-m)  condition  (this  is  in  the  region  north  of  Redfish 
Reef  and  along  Atkinson  Island).  In  this  upper  reach,  the  sedimentation 
dredging  volume  decline  is  moderately  steady  throughout  the  34  years  of 
this  channel  configuration.  In  the  lower  bay,  there  was  no  strong  trend 
over  this  period.  Analysis  of  the  freshwater  flow  to  the  system  for  the 
Trinity  and  San  Jacinto  rivers  appears  to  show  a  slight  increase  in  flow 
over  this  time  period.  However  there  is  not  a  strong  relationship  shown 
between  the  long-term  decline  in  required  dredging  and  the  freshwater  in¬ 
flow  to  the  system. 

It  is  possible  that  the  channel  stabilized  over  time  causing  less  sediment  to 
erode  from  its  banks.  Perhaps  dredging  practices  improved  so  that  less 
dredged  material  was  allowed  to  reenter  the  channel.  Another  possibility 
is  that  the  ship  traffic  has  increased  substantially  over  this  time  period. 
This  increase  would  have  an  impact  on  the  dredging  since  the  vessels  tend 
to  cause  resuspension  of  the  bed  material  meaning  that  less  sediment 
would  remain  on  the  bed.  Also,  a  possibility  for  the  drop  in  dredging  over 
this  period  is  the  change  in  the  sediment  load  from  upstream  sources  to 
the  bay.  Ward  and  Armstrong  (1993)  found  that  the  suspended  solids  in 
the  upper  Houston  Ship  Channel  had  dropped  to  one  third  of  that  from  25 
years  earlier.  They  attribute  this  reduction  to  improved  waste  treatment, 
altered  land  use,  and  impoundments  on  the  principal  rivers.  Solis  et.  al. 
(1994)  found  that  the  Trinity  River  gauge  at  Romayor,  TX  (USGS  gauge 
08066500  located  approximately  46  miles  (74  km)  upstream  from  the 
river’s  entrance  to  the  bay)  shows  a  significant  reduction  in  downstream 
sediment  inputs  after  the  completion  of  Lake  Livingston  in  1968.  There 
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actually  is  a  decline  in  total  suspended  sediment  throughout  most  of  the 
bay. 

Field  Data 

The  sediment  field  data  collection  revealed  finer  sediments  in  the  bed  of 
the  channel  and  at  a  considerable  distance  to  the  west,  and  also  to  the  east 
beyond  Atkinson  Island.  In  the  shallows  near  the  channel  the  bed  material 
was  firmer,  containing  more  sands.  The  area  east  of  the  channel  (toward 
Atkinson  Island)  contained  coarser  material  than  to  the  west,  which  is 
likely  a  result  of  exposure  to  the  waves  caused  by  vessel  traffic.  The  west¬ 
ern  edge  of  Atkinson  Island  contained  sand  and  shells.  Moving  around  the 
island  to  areas  more  protected,  the  bed  material  returned  to  finer  sedi¬ 
ments.  What  had  previously  been  islands  near  the  barge  mooring  basin 
had  eroded  to  about  a  6-ft  (i.8-m)  depth,  which  was  also  the  approximate 
depth  of  the  former  disposal  area  just  south  of  Atkinson  Island.  The  bed 
material,  which  was  quite  firm,  contained  progressively  more  fines  as  one 
moved  to  the  west  of  the  channel  out  into  the  bay.  This  distance  was 
greater  than  the  effect  of  the  vessels.  It  was  consistent  with  the  natural 
deposition  pattern  of  a  river  flowing  into  a  bay.  The  currents  slowed  down 
and  the  coarser  materials  settled  first.  Within  the  channel  (in  the  region 
around  Atkinson  Island)  the  bed  material  became  finer  moving  upland. 
The  currents  were  small  and  the  net  drift  near  the  bed  was  shown  to  be  up¬ 
land.  So  this  sediment  gradation  was  consistent  with  the  net  drift  in  the 
channel.  In  general,  the  suspended  sediment  concentrations  were  low 
during  the  survey.  Retrieving  samples  behind  vessels  did  not  necessarily 
result  in  high  concentrations.  The  highest  measurements  occurred  behind 
a  vessel  that  was  toward  the  edge  of  the  channel. 

Vessel  Effects 

A  vessel  moving  in  the  channel  causes  a  depression  in  the  water  surface 
and  return  currents  moving  from  the  bow  to  the  stern.  The  wavelength  of 
this  depression  is  roughly  that  of  the  vessel  length.  Generally  this  is  much 
greater  than  the  depth  of  the  waterway.  In  this  case,  the  wave  is  consid¬ 
ered  a  long  wave  and  travels  at  a  speed  of  about  ( ghf/ 2 .  This  speed  is 

called  the  wave  celerity.  This  depression  wave  moves  away  from  the  ves¬ 
sel.  The  speed  of  the  wave  is  calculated  using  the  local  depth.  If  the  vessel 
speed  is  less  than  the  wave  celerity,  then  the  vessel  impact  may  be  felt 
ahead  of  the  arrival  of  the  vessel.  Although  it  will  be  a  diminished  effect. 

If  the  vessel  is  traveling  in  the  channel  more  quickly  than  the  wave  speed 
in  the  shallows,  then  the  drawdown  will  arrive  after  the  vessel.  In  this  case 
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a  bore  can  be  formed.  With  a  channel  that  is  parallel  to  the  shoreline  (as  is 
the  case  with  the  Houston  Ship  Channel)  this  drawdown  is  reflected  at  the 
shoreline  and  by  the  deeper  water  forming  a  series  of  bores  perpendicular 
to  the  shore.  The  breaking  wave,  or  bore,  can  cause  high  shear  stresses 
and  erosion. 

This  effect  is  apparent  in  the  field  sediment  data  that  shows  the  grain  sizes 
on  the  Atkinson  Island  side  of  the  channel  to  be  larger  than  on  the  oppo¬ 
site  channel  side.  (The  shoreline  in  places  has  been  artificially  armored  as 
well.)  This  bed  material  is  now  coarser  than  the  usual  material  found  in 
channel  deposits.  The  increased  vessel  speed  associated  with  the  channel 
enlargement  does  contribute  to  increased  channel  deposition.  However, 
there  have  been  problems  with  erosion  of  the  islands  along  the  Houston 
Ship  Channel  for  some  time  before  the  channel  enlargement. 

A  vessel  moving  in  the  channel  will  cause  significant  return  currents  that 
can  cause  erosion  of  the  freshly  deposited  sediment.  It  was  estimated  that 
the  return  currents  should  decrease  in  the  enlarged  channel  even  though 
the  vessel  speed  increases.  Even  so,  the  amount  of  erosion  due  to  the  ves¬ 
sels  over  a  year  was  small  compared  to  the  annual  deposition.  It  does  not 
appear  that  at  this  point  the  decrease  in  direct  erosion  due  to  the  vessel 
could  be  responsible  for  a  large  increase  in  deposition  as  noted  in  the 
dredging  requirements. 

Hydrodynamic  and  Sediment  Model 

Generally  the  Houston  Ship  Channel  within  Galveston  Bay  is  flood- 
directed  bed  shear.  However,  for  the  base  conditions  in  the  months  of 
June  and  July  the  shear  is  ebb-directed  in  the  northern  portion  of  the  bay 
channel.  This  period  of  ebb-dominated  shear  is  likely  due  to  the  drop  in 
salinity  in  the  ocean  just  before  this  time.  The  dense  water  is  then  tending 
to  evacuate  the  system.  The  period  of  low  salinity  at  the  Gulf  of  Mexico  is 
due  to  the  freshwater  currents  from  the  Atchafalaya  and  Mississippi  Riv¬ 
ers.  For  the  other  months  and  for  the  plan  channel  the  bed  shear  is  flood- 
directed.  In  June  and  July  the  plan  condition  shears  are  flood-directed, 
but  much  less  strongly.  So  the  effect  of  the  drop  in  Gulf  salinity  still  re¬ 
mains,  but  the  larger  channel  prevents  the  complete  reversal  of  direction 
of  the  bed  shears. 

The  area  around  the  junction  of  the  Bayport  Channel  and  the  Houston 
Ship  Channel  shows  a  disrupted  shear  pattern.  The  shear  is  low  in  the 
flare  for  the  Bayport  Channel  but  also  in  the  Houston  Ship  Channel  near 
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that  location.  This  does  not  appear  when  the  shear  is  ebb-directed  but 
does  whenever  the  shear  is  flood-directed.  There  is  a  bend  in  the  Houston 
Ship  Channel  just  south  of  this  junction  that  may  be  the  cause. 

The  shear  stresses  throughout  the  upper  bay  are  slightly  lower  in  the  plan 
condition.  In  the  channel  the  plan  condition  showed  a  more  significant 
reduction  in  shear.  This  reduced  velocity  and  shear  will  result  in  a  greater 
deposition  potential. 

The  TABS-MDS  sediment  model  was  used  as  an  analysis  tool  to  determine 
where  suspended  sediment  travels  once  it  enters  the  channel  and  how  this 
pattern  changed  with  the  enlarged  channel.  A  single  grain  was  modeled 
and  an  equal  amount  of  sediment  was  introduced  into  both  channel  con¬ 
figurations  over  a  3-day  time  period.  This  simulation  was  conducted  for 
April,  May,  and  June  flow  conditions  -  each  month  modeled  separately. 
The  sediment  parameters  were  determined  from  the  field  data  collection. 

The  simulations  show  that  the  increase  in  sedimentation  in  the  channel  for 
the  larger  45-  x  530-ft  (13.7-  x  162-m)  channel  is  in  the  range  of  20- 
30  percent.  In  May  and  June  the  increase  is  concentrated  in  the  upland 
area  around  Atkinson  Island.  In  April  the  increase  is  concentrated  farther 
oceanward  in  the  Red  Fish  Reef  region.  This  simulation  does  not  include 
multiple  grains  or  sediment  entering  the  system  from  the  Trinity  or  San 
Jacinto  rivers. 

Conclusion 

This  investigation  considered  the  effects  of  vessels  on  erosion  in  and 
around  the  channel  and  the  change  in  currents  due  to  the  channel 
enlargement.  These  processes  will  change  the  channel  shoaling  long-term. 
Other  processes  such  as  leakage  from  disposal  sites,  channel  side  slope  re¬ 
laxation,  and  material  escape  during  dredging  affect  the  required  mainte¬ 
nance  dredging  over  a  short  time  period  (a  few  dredging  cycles).  Using  a 
sediment  tracer  and  results  from  field  data,  this  investigation  of  the  long¬ 
term  processes  only  revealed  a  potential  increase  in  dredging  on  the  order 
of  20  to  30  percent,  which  is  considerably  less  than  the  100  or  more  per¬ 
cent  increase  reported.  This  would  imply  that  the  large  increase  is  due  to  a 
short-term  process  and  is  not  likely  to  be  a  permanent  state.  These  results 
rely  upon  a  sediment  model  that  is  not  validated  to  match  field  results. 
Sediment  is  introduced  in  an  artificial  manner  as  a  tracer  and  not  through 
loads  in  the  river  system.  In  the  following  study  (Phase  2),  the  sediment 
model  will  be  validated  and  used  to  refine  the  predictions. 
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Appendix  A:  Field  Data  Analyses 

Tables  A1-A3  give  the  date,  time,  point  name,  and  State  Plane  coordinates 
for  all  of  the  data  points  taken  by  the  field  team  during  the  data  collection 
on  20-22  August  2004.  Following  these  tables  are  figures  showing  refer¬ 
ence  locations  for  all  of  the  data  points  (Figures  A1-A4).  Table  A4  gives 
the  suspended  concentration  and  salinity  for  many  of  the  field  data  points. 

Sediment  bulk  density  was  determined  using  a  pycnometer,  following 
standard  procedures.  The  results  of  this  analysis  are  given  in  Table  A5. 

Table  A6  along  with  Figures  A5  and  A6  show  results  of  the  Vertical  Loop 
Sediment-Water  Tunnel  (VOST)  analysis  and  Total  Suspended  Mass 
analysis.  A  detailed  explanation  of  the  VOST  can  be  found  in  Teeter  et  al. 
(1999).  The  VOST  tests  allow  for  the  determination  of  erosion  rate  for  the 
sediment  being  analyzed. 

Tables  A 7  through  A18  as  well  as  Figures  Ay  through  A16  present  results  of 
sediment  settling  rates.  These  tests  were  performed  in  the  standard 
method  of  determining  the  weight  of  sediment  on  the  bed  over  time  in  or¬ 
der  to  deduce  the  remaining  concentration  in  the  water  column. 

Standard  sieve  analyses  were  performed  on  push  core  #1  and  push  core 
#3.  These  results  are  presented  in  Tables  A19  and  A20. 

Figures  A17  through  A25  show  the  results  of  the  COULTER®  LS  Particle 
Size  Analyzer  for  several  field  data  points.  The  generated  plot  gives  a  vis¬ 
ual  image  of  the  grain  size  distribution  at  the  given  location. 

The  sand-silt/clay  separation  analysis  is  shown  in  Table  A21  and  Fig¬ 
ure  A26.  The  percentage  of  each  material  type  is  given  for  various  data 
points  around  the  channel  and  into  the  western  shallows. 
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Table  Al.  Field  data  point  information  for  20  August  2004. 


State  Plane  Coordinates 

NAD-83  Texas  South  Central  (4204),  US  survey  feet 

Date 

Point  Name 

X  coordinate 

Y  coordinate 

8/20/2004 

8:12:59 

3330723.27 

13847884.8 

8/20/2004 

9:00:09 

3242272.51 

13820918.2 

8/20/2004 

9:30:17 

3242050.23 

13820940.4 

8/20/2004 

9:30:18 

3242047.47 

13820943.1 

8/20/2004 

9:38:38 

3241411.21 

13820748.1 

8/20/2004 

9:43:19 

3241327.98 

13820640.6 

8/20/2004 

9:49:46 

3241329.88 

13820594.9 

8/20/2004 

9:51:32 

3241285.22 

13820565 

8/20/2004 

9:54:36 

3241389.03 

13820400.6 

8/20/2004 

10:03:41 

3242393.04 

13818555.1 

8/20/2004 

10:10:43 

3242913.22 

13819137.1 

8/20/2004 

10:13:45 

3242716.2 

13819388.8 

8/20/2004 

10:22:22 

3244030.36 

13818857.3 

8/20/2004 

10:29:55 

3245281.24 

13819169.3 

8/20/2004 

10:58:53 

3244320.15 

13815818.1 

8/20/2004 

11:03:27 

3244583.78 

13815865.2 

8/20/2004 

11:08:20 

3245043.95 

13816109.6 

8/20/2004 

11:10:42 

3245369.17 

13816356.4 

8/20/2004 

11:24:45 

3247491.07 

13811357.8 

8/20/2004 

11:36:23 

3247077.98 

13810943.9 

8/20/2004 

11:42:05 

3246419.76 

13810953.3 

8/20/2004 

11:45:54 

3245862.1 

13810894.6 

8/20/2004 

11:58:19 

3248906.4 

13807185.7 

8/20/2004 

12:04:47 

3248288.74 

13807290.6 

8/20/2004 

12:11:02 

3247773.12 

13807193.6 

8/20/2004 

12:14:18 

3247449.68 

13807060.3 

8/20/2004 

12:30:28 

3248867.81 

13802934.7 

8/20/2004 

12:34:17 

3249297.42 

13803056.7 

8/20/2004 

12:39:24 

3250113.36 

13803316.9 

8/20/2004 

12:43:01 

3249763.18 

13803177.4 

8/20/2004 

12:58:28 

3252139.99 

13798414.9 

8/20/2004 

13:01:34 

3251577.37 

13798328.3 

8/20/2004 

13:07:19 

3250963.1 

13798304.4 

8/20/2004 

13:07:21 

3250960.18 

13798304.3 

8/20/2004 

13:14:42 

3250862.83 

13797831.2 

8/20/2004 

13:38:06 

3255695.51 

13791238.4 
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8/20/2004 

13:43:04 

3254581.73 

13790652 

8/20/2004 

13:52:57 

3254134.45 

13790181.1 

8/20/2004 

13:58:30 

3253879.99 

13790108.6 

8/20/2004 

14:14:10 

3246023.46 

13804228.7 

8/20/2004 

14:36:42 

3248073.99 

13820713.5 

8/20/2004 

14:43:35 

3248213.44 

13820088.1 

8/20/2004 

14:47:59 

3248031.28 

13819930.9 

Table  A2.  Field  data  point  information  for  21  August  2004. 


State  Plane  Coordinates 

NAD-83  Texas  South  Central  (4204),  US  survey  feet 

Date 

Point  Name 

X  coordinate 

Y  coordinate 

8/21/2004 

10:50:11 

3242287.26 

13820978.4 

8/21/2004 

11:00:31 

3241888.14 

13820729.6 

8/21/2004 

11:13:54 

3242422.8 

13818552.9 

8/21/2004 

11:21:19 

3243508.28 

13818894.4 

8/21/2004 

11:26:07 

3244183.52 

13819120.3 

8/21/2004 

11:39:40 

3243414.15 

13819062.8 

8/21/2004 

11:48:04 

3244269.27 

13815853.6 

8/21/2004 

11:55:11 

3244588.61 

13816215.8 

8/21/2004 

12:00:34 

3245254.54 

13816409.4 

8/21/2004 

12:12:36 

3245418.22 

13813909.2 

8/21/2004 

12:19:47 

3246711.79 

13810646.7 

8/21/2004 

12:25:40 

3247103.04 

13810975.8 

8/21/2004 

12:40:33 

3247898.31 

13812336.6 

8/21/2004 

13:12:39 

3246677.68 

13810631 

8/21/2004 

13:18:05 

3246635.92 

13810796.4 

8/21/2004 

13:20:49 

3246122.84 

13810465.5 

8/21/2004 

13:27:01 

3246685.03 

13810614.6 

8/21/2004 

13:38:52 

3247735.36 

13810301.7 

8/21/2004 

13:46:00 

3248966.53 

13807163.4 

8/21/2004 

13:49:16 

3248284.92 

13807205.9 

8/21/2004 

13:52:20 

3247728.64 

13807185.8 

8/21/2004 

13:55:06 

3247426.76 

13807076.5 

8/21/2004 

14:00:16 

3248894.96 

13802916.9 

8/21/2004 

14:03:36 

3249289.17 

13803046.7 

8/21/2004 

14:08:07 

3248392.07 

13803763.6 

8/21/2004 

14:14:38 

3249808.53 

13803141.8 

8/21/2004 

14:19:10 

3250136.08 

13803368.5 
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8/21/2004 

14:31:56 

3252139.47 

13798449.7 

8/21/2004 

14:37:46 

3251514.64 

13798234 

8/21/2004 

14:45:15 

3251141.47 

13798334.5 

8/21/2004 

14:48:24 

3250837.6 

13798132.6 

8/21/2004 

15:06:52 

3251557.92 

13793036.5 

8/21/2004 

15:22:07 

3251551.74 

13793069.6 

8/21/2004 

15:33:38 

3254405.26 

13790458.9 

8/21/2004 

15:33:38 

3254405.26 

13790458.9 

8/21/2004 

15:44:34 

3259121.34 

13786438.8 

8/21/2004 

15:44:34 

3259121.34 

13786438.8 

8/21/2004 

16:09:26 

3262670.04 

13783473.6 

8/21/2004 

16:25:08 

3262696.94 

13780637.9 

8/21/2004 

16:30:07 

3262266.83 

13780324.8 

8/21/2004 

16:40:15 

3261840.1 

13780145.9 

Table  A3.  Field  data  point  information  for  22  August  2004. 


State  Plane  Coordinates 

NAD-83  Texas  South  Central  (4204),  US  survey  feet 

Date 

Point  Name 

X  coordinate 

Y  coordinate 

8/22/2004 

7:41:51 

3241917.58 

13802086 

8/22/2004 

8:06:00 

3243896.29 

13803519.9 

8/22/2004 

8:21:24 

3246521.38 

13805371.8 

8/22/2004 

8:48:36 

3257080.3 

13804196.1 

8/22/2004 

9:06:02 

3253510.73 

13807014.2 

8/22/2004 

9:26:35 

3255865.85 

13813601.1 

8/22/2004 

9:40:26 

3253121.33 

13818971.5 

8/22/2004 

10:02:35 

3248054.73 

13821831.8 

8/22/2004 

10:20:17 

3244940.34 

13819569.8 

8/22/2004 

11:02:25 

3238355.16 

13822505.4 

8/22/2004 

12:54:21 

3245970.58 

13815210.6 

8/22/2004 

13:08:07 

3248407.1 

13811501.9 

8/22/2004 

13:12:52 

3248986.2 

13807289.3 

8/22/2004 

13:18:46 

3250363.32 

13803790.6 

8/22/2004 

13:24:35 

3248624.06 

13802756.2 

8/22/2004 

13:29:21 

3247186.98 

13806829.2 

8/22/2004 

13:40:56 

3245774.1 

13810650.8 

8/22/2004 

13:52:43 

3244610.73 

13814099.2 
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Data  Points  for  20  August  2004 
(listed  West  to  East,  shown  on  45X530  ft  grid) 


Range  1  Points 

9:54:36 

9:51:32 

9:49:46 

9:43:19 

9:38:38 

9:30:18 

9:30:17 

9:00:09 


Range  2  Points 

10:03:41 

10:13:45 

10:10:43 

10:22:22 

10:29:55 

Range  3  Points 

10:58:53 

11:03:27 

11:08:20 

11:10:42 

Range  4  Points 

11:45:54 

11:42:05 

11:36:23 

11:24:45 

Range  5  Points 

12:14:18 

12:11:02 

12:04:47 

11:58:19 

Range  6  Points 

12:30:28 

12:34:17 

12:43:01 

12:39:24 


Range  7  Paints 

13.14:42 

13:07:21 

13:07:19 

13:01:34 

12:58:28 


Range  8  Points 

13:58:30 

13:52:57 

13:43:04 

13:38:06 


Other  Points 

#  14:14:10 

#  14:47  59 

#  14:36.42 

#  14:43:35 


Figure  Al.  Reference  map  for  20  August  2004  data  collection. 
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Data  Points  for  21  August  2004 
(listed  West  to  East,  shown  on  45X530  ft  grid) 


Range  1  Points 

11:00:31 

10:50:11 

Range  2  Points 

11:13:54 

11:39:40 

11:21:19 

11:26:07 


Range  3  Points 

11:48:04 

11:55:11 

12:00:34 

Range  4  Points 

13:20:49 

13:18:05 

13:12:39 

13:27:01 

12:19:47 

12:25:40 

13:38:52 

Range  5  Points 

13:55:06 

13:52:20 

13:49:16 

13:46:00 

Range  6  Points 

14:00:16 
14:03:36 
14  14:38 
14:19:10 


Range  7  Points 

14:48:24 

14:45:15 

14:37:46 

14.31:56 


Other  Points 

•12:12:36 
"12:40:33 
*"14:08:07 
tt  15:06:52 
#  152207 


Other  Points 

m  15:33:38  |  N 

$15:44:34 
$$16:09  26 
%  16:40:15 
%  16:30:07 
%  16:25:08 


Figure  A2.  Reference  map  for  21  August  2004  data  collection. 
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Data  Points  for  22  August  2004 
(shown  on  45X530  ft  grid) 


Figure  A3.  Reference  map  for  22  August  2004  data  collection. 
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Table  A4.  Suspended  sediment  sample  analysis  results. 


SUSPENDED  SEDIMENT  SAMPLE  ANALYSIS 

PROJECT:  Houston  Ship  Channel 

P.I.:  C.  Berger 

Sampling 

Range 

Sta.# 

Depth  Cone.  Salinity 

ft.  mg/I  ppt 

Date 

Time 

8/21/04 

105011 

1A 

6  O.B. 

19 

14.96 

8/21/04 

105011 

1A 

3  O.B. 

22 

15.03 

8/21/04 

110031 

1CL 

3  O.B. 

45 

15.43 

8/21/04 

110031 

1CL 

25  B.S. 

28 

15.15 

8/21/04 

110364 

2D 

3  B.S. 

18 

14.99 

8/21/04 

110364 

2D 

3  O.B. 

28 

14.99 

8/21/04 

112119 

2CL 

3  O.B. 

32 

15.57 

8/21/04 

112607 

2A 

3  O.B. 

24 

15.08 

8/21/04 

113940 

2CL 

3  O.B. 

49 

15.56 

8/21/04 

114804 

3D 

3  O.B. 

32 

15.18 

8/21/04 

115511 

3CL 

3  O.B. 

181 

15.56 

8/21/04 

120834 

3A 

3  B.S. 

53 

15.20 
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8/21/04 

121236 

3CL 

3  O.B. 

469 

15.90 

8/21/04 

121947 

4CL 

3  O.B. 

290 

16.42 

8/21/04 

121947 

4CL 

25  B.S. 

41 

15.87 

8/21/04 

122540 

3  O.B. 

25 

15.30 

8/21/04 

124033 

3  O.B. 

22 

15.22 

8/21/04 

131239 

3  O.B. 

107 

15.82 

8/21/04 

131805 

3  B.S. 

46 

15.67 

8/21/04 

132049 

3  O.B. 

35 

15.54 

8/21/04 

132704 

3  B.S. 

58 

15.65 

8/21/04 

132704 

3  O.B. 

2473 

16.17 

8/21/04 

132704 

Plume 

3  B.S. 

128 

15.67 

8/21/04 

133852 

3  O.B. 

32 

15.29 

8/21/04 

134600 

3  O.B. 

47 

15.44 

8/21/04 

134916 

3  O.B. 

145 

16.33 

8/21/04 

135220 

3  O.B. 

165 

16.23 

8/21/04 

135506 

3  O.B. 

44 

15.72 

8/21/04 

140006 

3  O.B. 

51 

15.89 

8/21/04 

140336 

3  O.B. 

93 

16.15 

8/21/04 

141438 

3  O.B. 

84 

16.19 

8/21/04 

141910 

3  O.B. 

46 

16.01 

8/21/04 

143156 

3  O.B. 

58 

15.96 

8/21/04 

143746 

3  O.B. 

97 

16.33 

8/21/04 

144515 

3  O.B. 

131 

16.34 

8/21/04 

144824 

3  O.B. 

45 

16.14 

8/21/04 

152207 

4  O.B. 

53 

17.84 

8/21/04 

153338 

3  O.B. 

138 

17.56 

8/21/04 

154434 

3  O.B. 

23 

13.61 

8/21/04 

160926 

4  B.S. 

67 

15.61 

8/21/04 

160926 

3  O.B. 

55 

14.99 

8/21/04 

162508 

3  B.S. 

46 

15.91 

8/21/04 

163007 

3  O.B. 

56 

16.25 

8/21/04 

164015 

3  O.B. 

202 

17.06 

8/22/04 

080600 

Tube  #2 

3  O.B. 

27 

15.64 

8/22/04 

080600 

Tube  #3 

3  O.B. 

73 

15.62 

8/22/04 

084836 

3  B.S. 

22 

11.52 

8/22/04 

092635 

3  O.B. 

29 

12.95 

8/22/04 

094026 

3  B.S. 

34 

14.72 

8/22/04 

100235 

2.5  B.S. 

32 

15.38 

8/22/04 

102017 

2.5  B.S. 

21 

15.31 

8/22/04 

110225 

3  B.S. 

33 

14.94 
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Table  A5.  Bulk  density  analysis  results  for  several  field  data  points. 


BULK  DENSITY  USING  PYCNO  BOTTLES 

PROJECT  NAME:  Houston  Ship  Channel 

Date  of  Project:  Aug.  04 

P.l. :  C.  Berger 

ANALYST:  D.B. 

Analysis  date: 

9/15/04 

Pycno. 

bottle  # 

Sample 

Date 

Sample 

ID 

Temp.  C 

Tare 

wt.  g 

Bottle  wt. 

+  sed.  g 

Bottle  wt. 

+  sed. 

+  water  g 

Density 

of 

water 

Volume 

Pycno. 

bottle 

Sediment 

Density 

g/cm3 

202 

8/20/04 

090009  10' 

22.5 

29.6822 

53.4343 

65.1580 

0.998 

24.9579 

1.798 

204 

8/20/04 

093838  48' 

22.5 

27.3518 

44.8146 

55.7142 

0.998 

25.6506 

1.186 

205 

8/20/04 

112445  6’ 

22.5 

28.3951 

50.5627 

64.0250 

0.998 

24.9066 

1.942 

208 

8/20/04 

114554  11' 

22.5 

27.3689 

49.8874 

61.9600 

0.998 

24.9452 

1.753 

209 

8/20/04 

121418  13' 

22.5 

27.9256 

50.4893 

63.2080 

0.998 

25.3391 

1.791 

202 

8/20/04 

123028  11' 

22.5 

29.6817 

50.1194 

63.0370 

0.998 

24.9579 

1.701 

204 

8/20/04 

123417  39’ 

22.5 

27.3517 

49.2058 

58.5816 

0.998 

25.6506 

1.344 

205 

8/20/04 

123924  9' 

22.5 

28.3950 

49.6617 

63.4850 

0.998 

24.9066 

1.924 

208 

8/20/04 

124301 39' 

22.5 

27.3688 

50.8687 

59.4858 

0.998 

24.9452 

1.441 

209 

8/21/04 

163007  13' 

22.5 

27.9258 

44.1559 

57.1981 

0.998 

25.3391 

1.323 

202 

8/21/04 

164015  51' 

22.5 

29.6820 

46.0553 

57.1184 

0.998 

24.9579 

1.180 

204 

8/22/04 

090602  4' 

22.5 

27.3518 

46.5548 

59.2823 

0.998 

25.6506 

1.489 

205 

10/26/04 

8/22/04 

080600  9.5' 

(core  tube  #2) 

22.5 

28.3948 

46.0233 

59.5205 

0.998 

24.9066 

1.549 

202 

11/3/04 

8/22/04 

074151 10.5' 

(core  tube  #1) 

23.0 

29.6821 

45.7397 

59.1796 

0.998 

24.9579 

1.397 

202 

8/22/04 

082124  9' 

(core  tube  #3) 

22.2 

29.6820 

50.8004 

63.0060 

0.998 

24.9579 

1.659 
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Table  A6.  Push  core  sample  analysis  results. 


Project  Name:  Houston  Ship  Channel 

P.I.:  C.  Berger 

Date  of  Project:  Aug.  20-22-04 

Date  of  Analysis: 

Vloop  ran  10/19/04 

Sample  I.D.: 

push  core  sample  #1  surface  slice 

TSM  ran  10/25/04 

8/22/04  074151 10.5  ft. 

Salinity 

-15  ppt 

Gross  wt. 

Shear 

Time 

Sample 

Tare  wt. 

(g) 

Net  wt. 

Cone. 

Bottle  # 

Voltage 

Stress 

(min.) 

Volume 

(g) 

filter  +  sed 

(g) 

(g/l) 

l 

3 

0.35 

10 

99 

0.0149 

0.0150 

0.0001 

0.0010 

2 

4 

0.50 

10 

100 

0.0152 

0.0153 

0.0001 

0.0010 

3 

5 

0.67 

10 

99 

0.0151 

0.0152 

0.0001 

0.0010 

4 

6 

0.83 

10 

99 

0.0151 

0.0153 

0.0002 

0.0020 

5 

7 

0.98 

10 

98 

0.0147 

0.0149 

0.0002 

0.0020 

6 

8 

1.14 

10 

100 

0.0151 

0.0153 

0.0002 

0.0020 

7 

8 

1.14 

15 

100 

0.0151 

0.0157 

0.0006 

0.0060 

8 

8 

1.14 

20 

100 

0.0148 

0.0157 

0.0009 

0.0090 

9 

8 

1.14 

25 

98 

0.0148 

0.0157 

0.0009 

0.0092 

10 

8 

1.14 

30 

99 

0.0149 

0.0158 

0.0009 

0.0091 

11 

8 

1.14 

35 

100 

0.0152 

0.0160 

0.0008 

0.0080 

12 

8 

1.14 

40 

100 

0.0149 

0.0159 

0.0010 

0.0100 

13 

8 

1.14 

45 

99 

0.0152 

0.0161 

0.0009 

0.0091 

14 

8 

1.14 

50 

99 

0.0150 

0.0159 

0.0009 

0.0091 

15 

9 

1.29 

5 

100 

0.0150 

0.0159 

0.0009 

0.0090 

16 

9 

1.29 

10 

100 

0.0150 

0.0161 

0.0011 

0.0110 

17 

9 

1.29 

15 

99 

0.0149 

0.0162 

0.0013 

0.0131 

Project  Name:  Houston  Ship  Channel 

P.I.:  C.  Berger 

Date  of  Project:  Aug.  20-22-04 

Date  of  Analysis: 

Vloop  ran  10/28/04 

Sample  I.D.: 

push  core  sample  #3  surface  slice 

TSM  ran  11/1/04 

8/22/04  082124  9.0  ft. 

Salinity 

-15  ppt 

Gross  wt. 

Shear 

Time 

Sample 

Tare  wt. 

(g) 

Net  wt. 

Cone. 

Bottle  # 

Voltage 

Stress 

(min.) 

Volume 

(g) 

filter  +  sed 

(g) 

(g/l) 

l 

3 

0.35 

10 

97 

0.0148 

0.0149 

0.0001 

0.0010 

2 

4 

0.50 

10 

99 

0.0145 

0.0146 

0.0001 

0.0010 

3 

5 

0.67 

10 

97 

0.0146 

0.0147 

0.0001 

0.0010 

4 

6 

0.83 

10 

96 

0.0148 

0.0149 

0.0001 

0.0010 

5 

7 

0.98 

10 

97 

0.0147 

0.0148 

0.0001 

0.0010 

6 

8 

1.14 

10 

100 

0.0147 

0.0148 

0.0001 

0.0010 

7 

9 

1.29 

10 

100 

0.0148 

0.0149 

0.0001 

0.0010 

8 

10 

1.45 

10 

97 

0.0148 

0.0149 

0.0001 

0.0010 
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9 

11 

1.62 

10 

98 

0.0148 

0.0151 

0.0003 

0.0031 

10 

12 

1.79 

10 

98 

0.0156 

0.0159 

0.0003 

0.0031 

11 

12 

1.79 

15 

100 

0.0156 

0.0158 

0.0002 

0.0020 

12 

13 

1.98 

5 

99 

0.0156 

0.0160 

0.0004 

0.0040 

13 

13 

1.98 

10 

98 

0.0155 

0.0161 

0.0006 

0.0061 

14 

13 

1.98 

15 

100 

0.0157 

0.0166 

0.0009 

0.0090 

15 

14 

2.15 

5 

100 

0.0154 

0.0172 

0.0018 

0.0180 

16 

14 

2.15 

10 

100 

0.0155 

0.0181 

0.0026 

0.0260 

17 

14 

2.15 

15 

98 

0.0156 

0.0197 

0.0041 

0.0418 

18 

14 

2.15 

20 

98 

0.0156 

0.0197 

0.0041 

0.0418 
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Push  Core  #1 


Volume  of  Water  (L)  =  6 

Surface  Area  of  Sample  (cm2)=  20.2683 


Cone. 

(9/1) 

Cone. 

(mg/I) 

Time 

(min.) 

Total  Time 
(min.) 

Shear 

Stress 

0.0010 

1.01 

10 

10 

0.35 

0.0010 

1.00 

10 

20 

0.50 

0.0010 

1.01 

10 

30 

0.67 

0.0020 

2.02 

10 

40 

0.83 

0.0020 

2.04 

10 

50 

0.98 

0.0020 

2.00 

10 

60 

1.14 

0.0060 

6.00 

15 

75 

1.14 

0.0090 

9.00 

20 

95 

1.14 

0.0092 

9.18 

25 

120 

1.14 

0.0091 

9.09 

30 

150 

1.14 

0  0080 

8.00 

35 

185 

1.14 

0.0100 

10.00 

40 

225 

1.14 

0.0091 

9.09 

45 

270 

1.14 

0.0091 

9.09 

50 

320 

1.14 

00090 

9  00 

5 

325 

1.29 

0  0110 

11.00 

10 

335 

1.29 

0.0131 

13.13 

15 

350 

1.29 

Shear 

Stress 

Change  in  Change  in 
Cone  Mass 

Erosion 

Erosion  Rate 

(Pa) 

(mg/L) 

(mg) 

Time  (mm) (mg/cm  mm) 

0.6700 

0.0100 

0.06 

30 

9  87E-05 

0.9800 

0.0900 

0.54 

30 

0.000888 

1.1400 

7  0500 

42  3 

260 

0.008027 

1.2900 

4.1300 

24.78 

25 

0.048904 

Ctfw  Ssmolti 


•  •  W 


Erosion  Rate  -  Push  Core  Sample  1 


0  2  0  4  0  6  0  8  1 

Shear  Stress  (Pa) 


1.2 


1.4 


Figure  A5.  Erosion  rate  analysis  for  push  core  #1. 
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Push  Core  #3 


Volume  of  Water  (L)  =  6 

Surface  Area  of  Sample  (cm2)=  20.2683 


Cone. 

Cone. 

Time 

Total  Time 

Shear 

O/l) 

(mg/I) 

(min.) 

(min.) 

Stress 

0.0010 

1.03 

10 

10 

0.35 

0.0010 

1.01 

10 

20 

0.50 

0.0010 

1.03 

10 

30 

0.67 

0.0010 

1.04 

10 

40 

0.83 

0.0010 

1.03 

10 

50 

0.98 

0.0010 

1.00 

10 

60 

1.14 

0.0010 

1.00 

10 

70 

1.29 

0.0010 

1.03 

10 

80 

1.45 

0.0031 

3.06 

10 

90 

1.62 

0.0031 

3.06 

10 

100 

1.79 

0.0020 

2.00 

15 

115 

1.79 

0  0040 

4.04 

5 

120 

1.98 

00061 

6.12 

10 

130 

1.98 

0.0090 

9.00 

15 

145 

1.98 

0  0180 

18.00 

5 

150 

2.15 

0.0260 

26.00 

10 

160 

2.15 

0  0418 

41  84 

15 

175 

2.15 

0.0418 

41.84 

20 

195 

2.15 

Shear 

Change  in 

Change  in 

Erosion 

Stress 

(Pa) 

Cone. 

(mg/L) 

Mass 

(mg) 

Time 

(min) 

(mg/cm2*min) 

1  7900 

2.0300 

12.18 

35 

0.01717 

1  9800 

5.9400 

35.64 

30 

0.058614 

2  1500 

32  8400 

197.04 

50 

0.194432 

Erosion  Rate  -  Push  Core  Sample  3 


Shear  Stress  (Pa) 


Figure  A6.  Erosion  rate  analysis  for  push  core  #3. 
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Table  A7.  Settling  test  results  for  field  data  point  114554,  run  1. 


SETTLING  TESTS-SUSPENDED  SEDIMENT  ANALYSIS 

1 

PROJECT:  Houston  Ship  Channel 

Instant  Ocean  water  ~15ppt 

Sample  I.D.:  114554  11  ft. 

12/8/04 

Cone.: 

3.4  ml/5  gal.  Cone.  target:50  mg/I 

Analyst:DB 

D.B 

Height  of  water  in  column=186.0  cm 

BWD=1.753 

Vol 

Weight,  grams 

Cone 

Time 

Notes 

ml 

Tare 

Gross 

Net 

mg/I 

0  min 

jerrican 

50 

0.0156 

0.0180 

0.0024 

48 

0  min 

column 

45 

0.0154 

0.0186 

0.0032 

71 

20 

51 

0.0156 

0.0190 

0.0034 

67 

40 

52 

0.0156 

0.0183 

0.0027 

52 

60 

51 

0.0157 

0.0180 

0.0023 

45 

90 

50 

0.0158 

0.0179 

0.0021 

42 

120 

51 

0.0156 

0.0175 

0.0019 

37 

150 

47 

0.0157 

0.0173 

0.0016 

34 

200 

52 

0.0156 

0.0173 

0.0017 

33 

250 

50 

0.0159 

0.0175 

0.0016 

32 

300 

46 

0.0159 

0.0173 

0.0014 

30 

350 

50 

0.0156 

0.0171 

0.0015 

30 
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Table  A8.  Settling  test  results  for  field  data  point  114554,  run  2. 


SETTLING  TESTS-SUSPENDED  SEDIMENT  ANALYSIS 

2 

PROJECT:  Houston  Ship  Channel 

Instant  Ocean  water  ~15ppt 

Sample  I.D.:  114554  lift. 

12/9/04 

Cone.: 

13.6  ml/5  gal.  Cone.  target:200 
mg/I 

Analyst: 

D.B. 

Height  of  water  in  column=186.0  cm 

BWD=1.753 

Vol 

Weight,  grams 

Cone 

Time 

Notes 

ml 

Tare 

Gross 

Net 

mg/I 

0  min 

jerrican 

50 

0.0155 

0.0248 

0.0093 

186 

0  min 

column 

47 

0.0155 

0.0247 

0.0092 

196 

20 

50 

0.0158 

0.0229 

0.0071 

142 

40 

50 

0.0157 

0.0217 

0.0060 

120 

60 

53 

0.0157 

0.0214 

0.0057 

108 

90 

46 

0.0155 

0.0201 

0.0046 

100 

120 

45 

0.0156 

0.0200 

0.0044 

98 

150 

50 

0.0155 

0.0203 

0.0048 

96 

200 

50 

0.0157 

0.0200 

0.0043 

86 

250 

50 

0.0158 

0.0196 

0.0038 

76 

300 

50 

0.0156 

0.0186 

0.0030 

60 

350 

47 

0.0158 

0.0181 

0.0023 

49 
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Table  A9.  Settling  test  results  for  field  data  point  114554,  run  3. 


SETTLING  TESTS-SUSPENDED  SEDIMENT  ANALYSIS 

3 

PROJECT:  Houston  Ship  Channel 

Instant  Ocean  water  ~15ppt 

Sample  I.D.:  114554  lift. 

12/13/04 

Cone.: 

13.6  +  20.6=  34.2  ml/5  gal.  Cone,  target:  500  mg/I  Analyst: 

D.B. 

Height  of  water  in  column=186.0  cm 

BWD=1.753 

Vol 

Weight,  grams 

Cone 

Time 

Notes 

ml 

Tare 

Gross 

Net 

mg/I 

0  min 

jerrican 

50 

0.0156 

0.0603 

0.0447 

894 

0  min 

column 

51 

0.0158 

0.0503 

0.0345 

676 

20 

51 

0.0157 

0.0339 

0.0182 

357 

40 

53 

0.0156 

0.0309 

0.0153 

289 

60 

50 

0.0156 

0.0291 

0.0135 

270 

90 

48 

0.0158 

0.0282 

0.0124 

258 

120 

47 

0.0156 

0.0264 

0.0108 

230 

150 

50 

0.0155 

0.0233 

0.0078 

156 

200 

50 

0.0157 

0.0208 

0.0051 

102 

250 

46 

0.0157 

0.0193 

0.0036 

78 

300 

50 

0.0159 

0.0189 

0.0030 

60 

350 

54 

0.0158 

0.0183 

0.0025 

46 
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Table  A10.  Settling  test  results  for  field  data  point  163007,  run  1. 


SETTLING  TESTS-SUSPENDED  SEDIMENT  ANALYSIS  1 

PROJECT:  Houston  Ship  Channel  Instant  Ocean  water  ~15ppt 

Sample  I.D.:  163007  13  ft.  12/8/04 

3.4  ml/5  gal.  Cone.  target:50  mg/I 

Cone.:  Analyst:  D.B. 

Height  of  water  in  column=186.0  cm  BWD=1.323 

Time 

Notes 

Vol 

ml 

Weight,  grams 

Cone 

mg/I 

Tare 

Gross 

Net 

0  min 

jerrican 

50 

0.0156 

0.0200 

0.0044 

88 

0  min 

column 

56 

0.0156 

0.0213 

0.0057 

102 

20 

53 

0.0155 

0.0217 

0.0062 

117 

40 

51 

0.0154 

0.0201 

0.0047 

92 

60 

54 

0.0155 

0.0200 

0.0045 

83 

90 

51 

0.0155 

0.0197 

0.0042 

82 

120 

56 

0.0156 

0.0199 

0.0043 

77 

150 

54 

0.0156 

0.0196 

0.0040 

74 

200 

56 

0.0155 

0.0196 

0.0041 

73 

250 

56 

0.0156 

0.0194 

0.0038 

68 

300 

52 

0.0156 

0.0188 

0.0032 

62 

350 

51 

0.0156 

0.0185 

0.0029 

57 
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Table  All.  Settling  test  results  for  field  data  point  163007,  run  2. 


SETTLING  TESTS-SUSPENDED  SEDIMENT  ANALYSIS 

2 

PROJECT:  Houston  Ship  Channel 

Instant  Ocean  water  ~15ppt 

Sample  I.D.:  163007  13  ft. 

12/9/04 

Cone.: 

13.6  ml/5  gal.  Cone.  target:200  mg/I 

Analyst: 

D.B. 

Height  of  water  in  column=186.0  cm 

BWD=1.323 

Vol 

Weight,  grams 

Cone 

Time 

Notes 

ml 

Tare 

Gross 

Net 

mg/I 

0  min 

jerrican 

50 

0.0155 

0.0253 

0.0098 

196 

0  min 

column 

48 

0.0155 

0.0280 

0.0125 

260 

20 

55 

0.0154 

0.0281 

0.0127 

231 

40 

50 

0.0154 

0.0263 

0.0109 

218 

60 

52 

0.0153 

0.0263 

0.0110 

212 

90 

49 

0.0155 

0.0258 

0.0103 

210 

120 

53 

0.0155 

0.0264 

0.0109 

206 

150 

52 

0.0154 

0.0245 

0.0091 

175 

200 

50 

0.0156 

0.0223 

0.0067 

134 

250 

51 

0.0156 

0.0202 

0.0046 

90 

300 

54 

0.0156 

0.0191 

0.0035 

65 

350 

54 

0.0156 

0.0185 

0.0029 

54 
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Table  A12.  Settling  test  results  for  field  data  point  163007,  run  3. 


SETTLING  TESTS-SUSPENDED  SEDIMENT  ANALYSIS 

3 

PROJECT:  Houston  Ship  Channel 

Instant  Ocean  water  15ppt 

Sample  I.D.:  163007  13  ft. 

12/13/04 

Cone.: 

13.6  +  20.6=  34.2  ml/5  gal.  Cone,  target:  500  mg/I  Analyst: 

D.B. 

Height  of  water  in  column=186.0  cm 

BWD=1.323 

Vol 

Weight,  grams 

Cone 

Time 

Notes 

ml 

Tare 

Gross 

Net 

mg/I 

0  min 

jerrican 

50 

0.0157 

0.0489 

0.0332 

664 

0  min 

column 

48 

0.0156 

0.0529 

0.0373 

777 

20 

50 

0.0156 

0.0546 

0.0390 

780 

40 

56 

0.0156 

0.0601 

0.0445 

795 

60 

53 

0.0156 

0.0576 

0.0420 

792 

90 

53 

0.0155 

0.0534 

0.0379 

715 

120 

55 

0.0158 

0.0398 

0.0240 

436 

150 

56 

0.0157 

0.0271 

0.0114 

204 

200 

54 

0.0156 

0.0217 

0.0061 

113 

250 

53 

0.0155 

0.0199 

0.0044 

83 

300 

55 

0.0158 

0.0195 

0.0037 

67 

350 

54 

0.0154 

0.0185 

0.0031 

57 
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Table  A13.  Settling  test  results  for  field  data  point  090009,  run  1. 


SETTLING  TESTS-SUSPENDED  SEDIMENT  ANALYSIS 

1 

PROJECT:  Houston  Ship  Channel 

Instant  Ocean  water  ~15ppt 

Sample  I.D.:  090009  10  ft. 

12/14/04 

Cone.: 

3.4  ml/5  gal.  Cone.  target:50  mg/I  Analyst: 

D.B. 

Height  of  water  in  column=186.0  cm 

BWD=1.798 

Vol 

Weight,  grams 

Cone 

Time 

Notes 

ml 

Tare 

Gross 

Net 

mg/I 

0  min 

jerrican 

50 

0.0156 

0.0180 

0.0024 

48 

0  min 

column 

48 

0.0156 

0.0193 

0.0037 

77 

20 

50 

0.0156 

0.0169 

0.0013 

26 

40 

53 

0.0156 

0.0168 

0.0012 

23 

60 

48 

0.0156 

0.0166 

0.0010 

21 

90 

48 

0.0156 

0.0166 

0.0010 

21 

120 

47 

0.0157 

0.0166 

0.0009 

19 

150 

48 

0.0156 

0.0165 

0.0009 

19 

200 

47 

0.0157 

0.0166 

0.0009 

19 

250 

46 

0.0155 

0.0163 

0.0008 

17 

300 

51 

0.0155 

0.0163 

0.0008 

16 

350 

50 

0.0156 

0.0164 

0.0008 

16 
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Table  A14.  Settling  test  results  for  field  data  point  090009,  run  2. 


SETTLING  TESTS-SUSPENDED  SEDIMENT  ANALYSIS 

2 

PROJECT:  Houston  Ship  Channel 

Instant  Ocean  water 

~15ppt 

Sample  I.D.:  090009  10  ft. 

12/15/04 

Cone.: 

13.6  ml/5  gal.  Cone.  target:200 
mg/I 

Analyst: 

D.B. 

Height  of  water  in  column=186.0  cm 

BWD=1.798 

Vol 

Weight,  grams 

Cone 

Time 

Notes 

ml 

Tare 

Gross 

Net 

mg/I 

0  min 

jerrican 

50 

0.0155 

0.0217 

0.0062 

124 

0  min 

column 

47 

0.0157 

0.0212 

0.0055 

117 

20 

50 

0.0153 

0.0195 

0.0042 

84 

40 

52 

0.0155 

0.0194 

0.0039 

75 

60 

48 

0.0156 

0.0187 

0.0031 

65 

90 

47 

0.0154 

0.0182 

0.0028 

60 

120 

46 

0.0158 

0.0185 

0.0027 

59 

150 

48 

0.0159 

0.0186 

0.0027 

56 

200 

49 

0.0154 

0.0181 

0.0027 

55 

250 

49 

0.0156 

0.0181 

0.0025 

51 

300 

46 

0.0158 

0.0181 

0.0023 

50 

350 

51 

0.016 

0.0184 

0.0024 

47 
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Table  A15.  Settling  test  results  for  field  data  point  090009,  run  3. 


SETTLING  TESTS-SUSPENDED  SEDIMENT  ANALYSIS 

3 

PROJECT:  Houston  Ship  Channel 

Instant  Ocean  water  ~15ppt 

Sample  I.D.:  090009  10  ft. 

12/16/04 

Cone.: 

13.6  +  20.6=  34.2  ml/5  gal.  Cone,  target:  500  mg/I  Analyst: 

D.B. 

Height  of  water  in  column=186.0  cm 

BWD=1.798 

Vol 

Weight,  grams 

Cone 

Time 

Notes 

ml 

Tare 

Gross 

Net 

mg/I 

0  min 

jerrican 

50 

0.0154 

0.0306 

0.0152 

304 

0  min 

column 

47 

0.0152 

0.0261 

0.0109 

232 

20 

51 

0.0154 

0.0247 

0.0093 

182 

40 

47 

0.0151 

0.0225 

0.0074 

157 

60 

51 

0.0155 

0.0228 

0.0073 

143 

90 

51 

0.0153 

0.0223 

0.0070 

137 

120 

46 

0.0155 

0.0216 

0.0061 

133 

150 

48 

0.0155 

0.0216 

0.0061 

127 

200 

50 

0.0153 

0.0208 

0.0055 

110 

250 

50 

0.0154 

0.0194 

0.0040 

80 

300 

47 

0.0156 

0.0185 

0.0029 

62 

350 

49 

0.0153 

0.0180 

0.0027 

55 
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Table  A16.  Settling  test  results  for  field  data  point  090602,  run  1. 


SETTLING  TESTS-SUSPENDED  SEDIMENT  ANALYSIS 

1 

PROJECT:  Houston  Ship  Channel 

Instant  Ocean  water  ~15ppt 

Sample  I.D.:  090602  4.0  ft. 

12/14/04 

Cone.: 

3.4  ml/5  gal.  Cone.  target:50  mg/I 

Analyst: 

D.B. 

Height  of  water  in  column=186.0  cm 

BWD=1.489 

Vol 

Weight,  grams 

Cone 

Time 

Notes 

ml 

Tare 

Gross 

Net 

mg/I 

0  min 

jerrican 

50 

0.0156 

0.0209 

0.0053 

106 

0  min 

column 

48 

0.0157 

0.0235 

0.0078 

163 

20 

52 

0.0156 

0.0219 

0.0063 

121 

40 

54 

0.0156 

0.0194 

0.0038 

70 

60 

53 

0.0155 

0.0184 

0.0029 

55 

90 

49 

0.0157 

0.0183 

0.0026 

53 

120 

50 

0.0157 

0.0181 

0.0024 

48 

150 

54 

0.0158 

0.0183 

0.0025 

46 

200 

53 

0.0154 

0.0178 

0.0024 

45 

250 

57 

0.0155 

0.0180 

0.0025 

44 

300 

52 

0.0153 

0.0177 

0.0024 

46 

350 

60 

0.0154 

0.0182 

0.0028 

47 
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Table  A17.  Settling  test  results  for  field  data  point  090602,  run  2. 


SETTLING  TESTS-SUSPENDED  SEDIMENT  ANALYSIS  2 

Instant  Ocean  water 

PROJECT:  Houston  Ship  Channel  ~15ppt 

Sample  I.D.:  090602  4.0  ft.  12/15/04 

13.6  ml/5  gal.  Cone.  target:200  mg/I 

Cone.:  Analyst:  D.B. 

Height  of  water  in  column=185.0  cm  BWD=1.489 

Time 

Notes 

Vol 

ml 

Weight,  grams 

Cone 

mg/I 

Tare 

Gross 

Net 

0  min 

jerrican 

50 

0.0154 

0.0334 

0.0180 

360 

0  min 

column 

52 

0.0157 

0.0389 

0.0232 

446 

20 

50 

0.0155 

0.0363 

0.0208 

416 

40 

59 

0.0159 

0.0324 

0.0165 

280 

60 

57 

0.0152 

0.0284 

0.0132 

232 

90 

51 

0.0154 

0.0267 

0.0113 

222 

120 

53 

0.0155 

0.0266 

0.0111 

209 

150 

53 

0.0159 

0.0257 

0.0098 

185 

200 

50 

0.0158 

0.0232 

0.0074 

148 

250 

52 

0.0152 

0.0205 

0.0053 

102 

300 

51 

0.0155 

0.0193 

0.0038 

75 

350 

59 

0.0155 

0.0191 

0.0036 

61 
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Table  A18.  Settling  test  results  for  field  data  point  090602,  run  3. 


SETTLING  TESTS-SUSPENDED  SEDIMENT  ANALYSIS 

3 

PROJECT:  Houston  Ship  Channel 

Instant  Ocean  water  ~15ppt 

Sample  I.D.:  090602  4.0  ft. 

12/16/04 

Cone.: 

13.6  +  20.6=  34.2  ml/5  gal.  Cone,  target:  500  mg/I  Analyst: 

D.B. 

Height  of  water  in  column=186.0  cm 

BWD=1.489 

Vol 

Weight,  grams 

Cone 

Time 

Notes 

ml 

Tare 

Gross 

Net 

mg/I 

0  min 

jerrican 

51 

0.0154 

0.0737 

0.0583 

1143 

0  min 

column 

48 

0.0153 

0.0649 

0.0496 

1033 

20 

53 

0.0154 

0.0597 

0.0443 

836 

40 

55 

0.0155 

0.0521 

0.0366 

665 

60 

52 

0.0154 

0.0449 

0.0295 

567 

90 

51 

0.0153 

0.0386 

0.0233 

457 

120 

53 

0.0153 

0.0304 

0.0151 

285 

150 

52 

0.0156 

0.0242 

0.0086 

165 

200 

56 

0.0153 

0.0218 

0.0065 

116 

250 

52 

0.0155 

0.0199 

0.0044 

85 

300 

55 

0.0157 

0.0192 

0.0035 

64 

350 

54 

0.0152 

0.0180 

0.0028 

52 

Table  18.  Settling  test  results  for  field  data  point  090602,  run  3. 
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090009  1 


Settling  Chamber  Experiment  Results  and  Settling  Velocity  Determination 

20-Jun-Q5 
Jennifer  N  Tato 


INPUT: 


CALCULATE: 


Time 

(mln) 

Concentration 

C0) 

Fraction  of 
Concentration 
Removed 

Ln  (Time) 

1 

0 

0048 

0 

2 

0 

0,077083333 

# 

n 

3 

20 

0.026 

0458333333 

2.9957323 

4 

40 

0  022641509 

0  528301887 

36888795 

5 

60 

0020833333 

0  565972222 

4  0943446 

6 

90 

0,020833333 

0.565972222 

4.4998097 

7 

120 

0,019148936 

0.60106383 

4  7874917 

8 

150 

001875 

0609375 

5.0106353 

9 

200 

0,019148936 

0  60106383 

5  2983174 

10 

250 

0  017391304 

0.637681 15S 

5  5214609 

11 

300 

0  015886775 

0  673202614 

5  7037825 

12 

350 

0016 

0666666667 

5  8579332 

13 

# 

* 

0 

n 

14 

# 

0 

# 

0 

15 

0 

0 

# 

0 

16 

# 

# 

0 

0 

17 

0 

# 

0 

18 

0 

# 

0 

0 

19 

0 

0 

0 

0 

20 

0 

0 

n 

0 

21 

# 

0 

n 

0 

22 

0 

0 

n 

0 

23 

* 

0 

n 

0 

24 

0 

# 

# 

0 

25 

0 

0 

0 

0 

26 

0 

0 

n 

# 

27 

0 

# 

# 

0 

28 

# 

# 

# 

0 

29 

# 

0 

« 

0 

30 

# 

0 

« 

0 

PuflHCar*  Siwplw 
intw# 


if  : 

/"l- 


Slope 

Y-intercept 

Roots 

0  0680755 

0  26768812 

-3  5102 

Chart  1 

A 

B 

C 

-0.03 

0  3417 

-0  3871 

mm/sec 

z-values 

b- 

dectile 

check 

Ws  Dectiles 

0.1 

0  05830689 

OK 

5.833 

0.2 

0.04630689 

OK 

3.764 

03 

0  03430689 

OK 

2.284 

04 

0  02230689 

OK 

1  256 

05 

0  01030689 

OK 

0  566 

06 

-0  0016931 

Imaginary 

0NUMI 

07 

0  0136931 

Imaginary 

0NUM1 

08 

-0  0256931 

Imaginary 

0NUM1 

09 

-0  0376931 

Imaginary 

0NUM1 

Chart  2 

- 7! - 

B 

C 

-0  0453 

0.4138 

•0  3757 

z-values 

b  = 

dectile 

check 

Ws  Dectiles 

6  1 

0  0850336 

OK 

8  047 

0.2 

0  0669136 

OK 

5  595 

0.3 

00487536 

OK 

3  687 

04 

0  0306736 

OK 

2  225 

05 

0  0125536 

OK 

1,109 

0.6 

-0.0055664 

Imaginary 

0NUMI 

0.7 

-0  0236864 

Imaginary 

0NUM1 

0.8 

-0.0418064 

Imaginary 

#NUM! 

09 

-0  0599264 

Imaginary 

0NUM! 

Figure  A7.  Settling  chamber  analysis  for  determination  of  settling  velocity  for  Point  090009, 

run  1. 
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090009  2 


Settling  Chamber  Experiment  Results  and  Settling  Velocity  Determination 


20-Jun-05 
Jennifer  N  Tato 


INPUT:  CALCULATE: 


Time 

(min) 

Concentration 

(9/1) 

Fraction  of 
Concentration 
Removed 

Ln  (Time) 

1 

ft 

0  124 

ft 

ft 

2 

# 

0,117021277 

ft 

ft 

3 

20 

0.084 

0  322580645 

29957323 

4 

40 

0  075 

0  39516129 

3  6888795 

5 

60 

0  064583333 

0  479166667 

4  0943446 

6 

90 

0,059574468 

0.519560741 

4.4998097 

7 

120 

0  058695652 

0  526647966 

4  7874917 

8 

150 

0.05625 

0.546370968 

5.0106353 

9 

200 

0,055102041 

0.555628703 

5  2983174 

10 

250 

0  051020408 

0.588545095 

5  5214609 

11 

300 

0  05 

0  596774194 

5  7037825 

12 

350 

0.047058824 

0  620493359 

5  8579332 

13 

ft 

ft 

ft 

ft 

14 

# 

# 

ft 

ft 

15 

# 

# 

ft 

ft 

16 

# 

ft 

0 

ft 

17 

# 

4 

ft 

ft 

18 

ft 

ft 

ft 

ft 

19 

ft 

Of 

ft 

ft 

20 

# 

ft 

ft 

ft 

21 

# 

# 

ft 

ft 

22 

ft 

# 

ft 

ft 

23 

ft 

ft 

ft 

ft 

24 

ft 

ft 

ft 

ft 

25 

ft 

ft 

ft 

ft 

26 

# 

ft 

ft 

ft 

27 

ft 

ft 

ft 

ft 

28 

# 

ft 

ft 

ft 

29 

# 

ft 

ft 

ft 

30 

# 

ft 

ft 

ft 

Slope 

Y-mtercept  Roots 

0  0986879 

0  04673616  -0.29192 

Chart  1 

A 

B  C 

-0  0134  0  2114 

2-values  b  — 

-0  2501 

dectile 

check 

mm/sec 

Ws  Dectiles 

0.1 

0  0259246 

OK 

4  729 

0.2 

0.0205646 

OK 

2.452 

03 

0  0152046 

OK 

1  158 

04 

0  0098446 

OK 

0472 

05 

0  0044846 

OK 

0  142 

06 

-0  0008754 

Imaginary 

ftNUM! 

07 

0  0062354 

Imaginary 

ftNUM' 

0.8 

-0  0115954 

Imaginary 

ftNUMl 

09 

-0  0169554 

Imaginary 

ftNUM' 

Chart  2 

- 7 - 

B 

C 

0  0132 

0.0372 

0  0901 

z-values 

b  = 

dectile 

check 

Ws  Dectiles 

0  1 

0  00190656 

OK 

24.267 

02 

0  00718656 

OK 

5  114 

0.3 

0  01246656 

OK 

1  847 

04 

0  01774656 

OK 

0.816 

05 

0  02302656 

OK 

0.405 

0.6 

0  02830656 

OK 

0.217 

07 

0  03358656 

OK 

0123 

0.8 

0  03886656 

OK 

0.072 

09 

0  04414656 

OK 

0  044 

Figure  A8.  Settling  chamber  analysis  for  determination  of  settling  velocity  for  Point  090009, 

run  2. 
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090009  3 


Settling  Chamber  Experiment  Results  and  Settling  Velocity  Determination 

20-Jun-Q5 
Jennifer  N  Tate 


INPUT: 


CALCULATE: 


Time 

(min) 

Concentration 

(3/1) 

Fraction  of 
Concentration 
Removed 

Ln  (Time) 

1 

0 

0  304 

0 

0 

2 

0 

0,231914594 

0 

0 

3 

20 

0.182352941 

0  400154799 

29957323 

4 

40 

0  157446809 

0482082867 

36888795 

5 

60 

0 143137255 

0  529153767 

4  0943446 

6 

90 

0.137254902 

0.548503612 

4.4998097 

7 

120 

0  132608696 

0  563787185 

4  7874917 

8 

150 

0.127083333 

0.581962719 

5.0106353 

9 

200 

0.11 

0  638157895 

5  2983174 

10 

250 

0.00 

0.736042105 

5  5214609 

11 

300 

0  061702128 

0  797032475 

5  7037825 

12 

350 

0.055102041 

0.818743287 

5  8579332 

13 

* 

0 

0 

0 

14 

# 

0 

0 

0 

15 

# 

0 

0 

0 

16 

# 

0 

0 

0 

17 

# 

0 

0 

0 

18 

# 

0 

0 

0 

19 

# 

0 

0 

0 

20 

0 

0 

0 

0 

21 

# 

0 

0 

0 

22 

0 

0 

0 

0 

23 

0 

0 

0 

0 

24 

0 

0 

0 

0 

25 

0 

0 

0 

0 

26 

0 

0 

0 

0 

27 

0 

0 

# 

0 

28 

0 

0 

0 

0 

29 

0 

0 

0 

0 

30 

# 

0 

0 

0 

tot*  Car* 


Slope 

Y-mtercept 

Roots 

0 1392844 

-0  0513794 

0  829392 

Chart  1 

A 

B 

C 

-0  0235 

02905 

-0  3338 

mm/sec 

z-values 

b- 

dectile 

check 

Ws  Dectiles 

0  1 

0  04361305 

OK 

5.455 

0.2 

0.03421305 

OK 

3.283 

03 

0  02481305 

OK 

1  831 

04 

0  01541305 

OK 

0  900 

05 

0  00601305 

OK 

0  334 

06 

-0  003387 

Imaginary 

0NUMI 

07 

0012787 

Imaginary 

0NUM" 

0.6 

-0  022187 

Imaginary 

0NUM1 

09 

-0  031587 

Imaginary 

0NUM' 

Chart  2 

- K - 

B 

C 

-0  0311 

0  3338 

-0  3214 
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b  = 
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check 

Ws  Dectiles 

0  1 
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OK 

7.189 

02 

0  04656028 

OK 

4  648 

0.3 

0  03412028 

OK 

2821 

04 

0  02168028 

OK 

1  544 

05 
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OK 
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Figure  A9.  Settling  chamber  analysis  for  determination  of  settling  velocity  for  Point  090009, 

run  3. 
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090602  1 


Settling  Chamber  Experiment  Results  and  Settling  Velocity  Determination 


20-Jun-Q5 
Jonnifor  N  Tato 


INPUT:  CALCULATE: 


Time 

(min) 

Concentration 

(fl/l) 

Fraction  of 
Concentration 
Removed 

Ln  (Time) 

1 

* 

0.106 

# 

4 

2 

# 

0.1625 

4 

4 

3 

20 

0.121153846 

-0  142960813 

2.9957323 

4 

40 

007037037 

0  336128581 

36888795 

5 

60 

0054716981 

0  483802065 

4  0943446 

6 

90 

0.053061224 

0  49942241 

4.4998097 

7 

120 

0.048 

0  547169811 

4  7874917 

8 

150 

0  046296296 

0  563242488 

5.0106353 

9 

200 

0.045283019 

0.572801709 

5298304 

10 

250 

0  043B59649 

0.586223725 

5  5214609 

11 

300 

0  046153846 

0  564586357 

5  7037825 

12 

350 

0  046666667 

0559748428 

5  8579332 

13 

4 

ft 

4 

n 

14 

4 

ft 

4 

4 

15 

4 

ft 

4 

» 

16 

# 

4 

4 

4 

17 

# 

4 

4 

4 

18 

# 

ft 

4 

4 

19 

ft 

4 

4 

4 

20 

ft 

ft 

4 

4 

21 

ft 

# 
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Figure  A10.  Settling  chamber  analysis  for  determination  of  settling  velocity  for  Point  090602, 

run  1. 
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Figure  All.  Settling  chamber  analysis  for  determination  of  settling  velocity  for  Point  090602, 

run  2. 
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Figure  A12.  Settling  chamber  analysis  for  determination  of  settling  velocity  for  Point  090602, 

run  3. 
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Figure  A13.  Settling  chamber  analysis  for  determination  of  settling  velocity  for  Point  114554, 

run  1. 
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Figure  A14.  Settling  chamber  analysis  for  determination  of  settling  velocity  for  Point  114554, 

run  2. 
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Figure  A15.  Settling  chamber  analysis  for  determination  of  settling  velocity  for  Point  114554, 

run  3. 
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Figure  A16.  Settling  chamber  analysis  for  determination  of  settling  velocity  for  Point  163007, 

run  1. 
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2.380 

No.8 

2.000 

No.  10 

1.414 

No.14 

1.190 

No.  16 

1.000 

No.  18 

840 

No.20 

.590 

No.30 

.500 

No.  35 

.0000  .000  .000  100.000 

.420 

No.40 

.0015  .208  .208  99.792 

.297 

No. 50 

250 

No. 60 

.210 

No. 70 

.0117  1.621  1.829  98.171 

149 

No.  100 

.0153  2.120  3.950  96.050 

.125 

No.  120 

105 

No. 140 

074 

No. 200 

.3409  47  242  51.192  48  808 

.063 

No. 230 

.2064  28,603  79.795  20.205 

PAN 

.1458  20.205  100.000  .000 

TOTAL  WEIGHT  IN 

GRAMS  .7216 

D95  = 

1467  mm 

D90  =  ,1362  mm  D85  =  ,1265  mm 

D80  = 

1175  mm 

D75  =  ,1091  mm  D70  =  .1013  mm 

D65  = 

0941  mm 

D60  =  ,0873  mm  D55  -  .081 1  mm 

D50  = 

0753  mm 

D45  =  0724  mm  D40  =  0704  mm 

D35  = 

0685  mm 

D30  =  0666  mm  D25  =  .0647  mm 
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Table  A20.  Sieve  analysis  results  for  push  core  #3. 


SIEVE  ANALYSIS 

PROJECT:  Houston  Ship  Channel 

SAMPLE  ID:  Core  Tube  #3  1-1/4  to  2  in-  below  sample  surface 

FILE  IS:  Core  3 

WEIGHT  PERCENT 

SIEVE  SIZE  RETAINED  PERCENT  RETAINED  FINER 

MILLIMETERS  OR  NUMBER  IN  GRAMS  PARTIAL  TOTAL  BY  WEIGHT 

76.200 

3  in 

50.800 

2  in 

38.100 

1-1/2  in 

25.400 

1  in 

19.050 

3/4  in 

12.700 

1/2  in 

9.525 

3/8  in 

6.727 

No.3 

4760 

No.4 

4.000 

No. 5 

3.360 

No. 6 

2.380 

No.8 

2.000 

No  10 

1.414 

No.14 

1.190 

No. 16 

1.000 

No. 18 

840 

No.20 

.590 

No.30 

500 

No. 35 

0000  .000  .000  100.000 

.420 

No.40 

0077  132  132  99  868 

.297 

No.  50 

.250 

No.60 

.210 

No.  70 

0684  1.170  1.302  98.698 

149 

No.  100 

.1041  1.781  3  084  96.916 

125 

No.  120 

.105 

No. 140 

.074 

No.200 

3.3412  57.175  60.259  39.741 

063 

No.230 

1.2523  21.430  81.688  18.312 

PAN 

1.0701  18.312  100.000  000 

TOTAL  WEIGHT  IN  GRAMS  5.8438 

D95  = 

1455  mm 

D90  =  ,1369  mm  D85=  1288  mm 

D80  = 

1211  mm 

D75  =  ,1139  mm  D70  =  1072  mm 

D65  = 

1008  mm 

D60  =  .0948  mm  D55  =  .0892  mm 

D50  = 

0839  mm 

D45  =  .0789  mm  D40  =  .0742  mm 

D35  = 

0714  mm 

D30  =  .0688  mm  D25  =  .0662  mm 

D20  = 

0638  mm 
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Coulter®  LS  Particle  Size  Analyzer 


File  name: 

093838 

Group  ID: 

Houston  Ship  Channel 

Sample  ID: 

093838  48  ft  8-20-04 

Location: 

Range  1,  channel 

Operator: 

DB 

Comments: 

oxi.  3  times 

Optical  model: 

Fraunhofer  rfz 

LS  100Q 

Fluid  Module 

Differential  Volume  (Average)  (2  S.D.) 


Particle  Diameter  (mm) 


Volume  Statistics  (Arithmetic)  093838 


Calculations  from  0  375  urn  to  948.3  urn 


Volume: 

100% 

Mean: 

9.679  pm 

S.D.: 

8.424  pm 

Median: 

6.384  pm 

C.V.: 

87.0% 

D(3,2). 

2.852  pm 

Skewness: 

0.658  Right  Skewed 

Mode: 

21.69  pm 

Kurtosis: 

-0.947  Pltykurtic 

%< 

10 

25 

50 

75 

90 

jim 

1.012 

2.257 

6.384 

17.28 

22.67 

Figure  A17.  Coulter  counter  analysis  for  Point  093838. 
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Coulter®)  LS  Particle  Size  Analyzer 


File  name. 

110820 

Group  ID: 

Houston  Ship  Channel 

Sample  ID: 

110820  39  ft.  8-20-04 

Location: 

Range  3,  channel,  Morgan's 

Point 

Operator: 

DB 

Comments: 

oxi.  3  times 

Optical  model: 

Fraunhofer.rfe 

LS  100Q 

Fluid  Module 

Differential  Volume  (Average)  (2  S.D.) 


Particle  Diameter  (mm) 


Volume  Statistics  (Arithmetic)  1 1 0820 
Calculations  from  0.375  pm  to  948.3  pm 


Volume: 

Mean: 

Median: 

D(3,2): 

Mode: 


100% 

1 1 .30  jim 
5.914  pm 
2.776  urn 
4.047  urn 


S.D.: 

C.V.: 

Skewness: 

Kurtosis: 


12.44  pm 
110% 

1 .445  Right  Skewed 
1.251  Leptokurtic 


%<  10  25  50  75  90 

pm  0.990  2.164  5.914  16.43  31.62 


Figure  A18.  Coulter  counter  analysis  for  Point  110820. 
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Coulter®  LS  Particle  Size  Analyzer 


File  name: 

114205 

Group  ID: 

Houston  Ship  Channel 

Sample  ID: 
Basin 

114205  47  ft  8-20-04 

Location: 

Range  4,  channel,  near  Mooring 

Operator 

DB 

Comments: 

oxi.  3  times 

Oplical  model: 

Fraunhofer. rfz 

LS  100Q 

Fluid  Module 

Differential  Volume  (Average)  (2  S.D.) 


Particle  Diameter  (mm) 


Volume  Statistics  (Arithmetic)  1 14205 
Calculations  from  0.375  pm  to  948.3  pm 


Volume: 

Mean: 

Median: 

D(3,2): 

Mode: 


100% 
8.410  pm 
4.714  pm 
2.503  pm 
4.877  pm 


S.D.: 

C.V.: 

Skewness: 

Kurtosis: 


8.922  pm 
106% 

1.471  Right  Skewed 
1.417  Leptokurtic 


%<  10  25  50  75  90 

pm  0.925  1.940  4.714  11.80  22.91 


Figure  A19.  Coulter  counter  analysis  for  Point  114205. 
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Coulter®  LS  Particle  Size  Analyzer 


File  name: 

120447 

Group  ID: 

Houston  Ship  Channel 

Sample  ID: 

120447  37  ft.  8-20-04 

Location: 

Range  5,  channel 

Operator: 

DB 

Comments: 

oxi.  3  times 

Optical  model 

Fraunhofer.rfz 

LS  100Q 

Fluid  Module 

Differential  Volume  (Average)  (2  S.D.) 


Particle  Diameter  (mm) 


Volume  Statistics  (Arithmetic)  120447 
Calculations  from  0.375  urn  to  948.3  urn 


Volume: 

Mean: 

Median: 

D(3,2): 

Mode: 


100% 
68.61  pm 
55.97  pm 
4.623  pm 
116.3  pm 


S.D.: 

C.V.: 

Skewness: 

Kurtosis: 


67.74  pm 
98.7% 

0.899  Right  Skewed 
0.429  Leptokurtic 


%<  10  25  50  75  90 

pm  1.386  4.311  55.97  115.3  147.8 


Figure  A20.  Coulter  counter  analysis  for  Point  120447. 
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Coulter®  LS  Particle  Size  Analyzer 


File  name: 

123417 

Group  ID: 

Houston  Ship  Channel 

Sample  ID: 

123417  39  ft.  8-20-04 

Location: 

Range  6.  channel 

Operator: 

DB 

Comments: 

oxi.  3  times 

Optical  model: 

Fraunhofer.rfz 

LS  100Q 

Fluid  Module 

Differential  Volume  (Average)  (2  S.D.) 


Particle  Diameter  (mm) 


Volume  Statistics  (Arithmetic)  123417 
Calculations  from  0  375  urn  to  948.3  pm 


Volume 

Mean: 

100% 

57.19  pm 

S.D.: 

50.91  pm 

Median: 

55.14  pm 

C.V.: 

89% 

D(3,2): 

5.072  pm 

Skewness: 

0.483  Right  Skewed 

Mode: 

96.49  pm 

Kurtosis: 

-0.801  Platykurtic 

%< 

10 

25  50 

75  90 

pm 

1.567 

5.135  56.14 

96.51  126.2 

Figure  A21.  Coulter  counter  analysis  for  Point  123417. 
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Coulter®  LS  Particle  Size  Analyzer 


File  name: 

123028 

Group  ID. 

Houston  Ship  Channel 

Sample  ID: 

123028  11  ft.  8-20-04 

Location: 

Range  6,  West  shallow 

Operator: 

D8 

Comments: 

oxi.  3  times 

Optical  model’ 

Fraunhofer.rfz 

LS  100Q 

Fluid  Module 

Differential  Volume  (Average)  (2  S.D.) 


Particle  Diameter  (mm) 


Volume  Statistics  (Arithmetic)  123028 


Calculations  from  0.375  um  to  948.3  am 


Volume: 

100% 

Mean: 

93.64  pm 

S.D.: 

51.64  pm 

Median: 

94.32  pm 

C.V.: 

55.1% 

D(3,2): 

10.83  pm 

Skewness: 

0.085  Right  Skewed 

Mode: 

105.9  pm 

Kurtosis: 

0.022  Leptokurtic 

%< 

10 

25 

50 

75 

90 

jim 

4.558 

65.62 

94.32 

122.6 

163.0 

Figure  A22.  Coulter  counter  analysis  for  Point  123028. 
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Coulter®  LS  Particle  Size  Analyzer 


Filename:  164015 

Group  ID: 

Houston  Ship  Channel 

Sample  ID:  164015  51  ft.  8-21-04 

Operator:  DB 

Comments:  oxi.  3  times 

Location: 

Channel,  near  marsh  site 

Optical  model:  Fraunhofer. rfz 

LS  100Q  Fluid  Module 

Differential  Volume  (Average)  (2  S.D.) 


Particle  Diameter  (mm) 


Volume  Statistics  (Arithmetic)  1 6401 5 
Calculations  from  0  375  urn  to  948.3  urn 


Volume: 

100% 

Mean: 

11.11  pm 

S.D.: 

14.18  pm 

Median: 

5.314  pm 

C.V.: 

128% 

D(3,2): 

2.657  pm 

Skewness: 

2.460  Right  Skewed 

Mode: 

4.443  pm 

Kurtosis: 

7.947  Leptokurtic 

%< 

10 

25 

50 

75 

90 

|im 

0.959 

2.054 

5.314 

14.88 

29.56 

Channel  near 
marsh  site 


V 


Figure  A23.  Coulter  counter  analysis  for  Point  164015. 
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Coulter®  LS  Particle  Size  Analyzer 


File  name: 

163007 

Group  ID: 

Houston  Ship  Channel 

Sample  ID: 

163007  13  ft.  8-21-04 

Location: 

East  shallows,  near  marsh  site 

Operator: 

D8 

Comments: 

oxi.  3  times 

Optical  model: 

Fraunhofer  rfz 

LS  100Q 

Fluid  Module 

Differential  Volume  (Average)  (2  S.D.) 


Particle  Diameter  (mm) 


Volume  Statistics  (Arithmetic)  163007 
Calculations  from  0  375  urn  to  948.3  urn 


Volume: 

100% 

Mean: 

14.15  pm 

S.D.: 

22.42  pm 

Median: 

5.654  pm 

C.V.: 

158% 

D(3,2): 

2.784  pm 

Skewness: 

3.072  Right  Skewed 

Mode: 

4.443  pm 

Kurtosis: 

10.23  Leptokurtic 

%< 

10 

25 

50 

75 

90 

|im 

1.001 

2.182 

5.654 

15.97 

34.14 

•*  i 

J-V 


\, 


East  shallows  . 
near  marsh  sile 


Figure  A24.  Coulter  counter  analysis  for  Point  163007. 
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Coulter®  LS  Particle  Size  Analyzer 


File  name: 
Sample  ID: 
Operator: 
Comments: 
Optical  model: 
LS  100Q 


PCTube2 

080600  9.5  ft.  8-22-04 
DB 

push  core  tube  #2,  oxl.  3  times 
Fraunhofer  rfz 
Fluid  Module 


Group  ID.  Houston  Ship  Channel 

Location:  Western  Lobe  shallows 


Differential  Volume  (Average)  (2  S.D.) 


Particle  Diameter  (mm) 


Volume  Statistics  (Arithmetic) 

PCTube  2 

Calculations  from  0.375  pm  to 

948.3  pm 

'<  \ 

\ 

m£m\ 

Volume: 

100% 

Mean: 

60.22  pm 

S.D.: 

74.50  pm 

Western  , 

Median: 

36.64  pm 

C.V.: 

124% 

1  dhe  shallow^  y  • 

D(3,2): 

4.587  pm 

Skewness: 

2.377  Right  Skewed 

V  .  \ 

Mode: 

87.90  pm 

Kurtosis: 

9.156  Leptokurtic 

^  A 

%<  10 

25  50 

75  90 

pm  1.383 

4.476  36.64 

89.37  146.5 

• 

** 

Figure  A25.  Coulter  counter  analysis  for  Point  0806009. 
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Figure  A26.  Grain  class  separation  percentages  for  several  field  data  points. 
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